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ABSTRACT 

The  aim  of  this  investigation  was  to  review  the  different  types  of  nearshore  current, 
turbidity,  and  morphology  evolution  models  currently  available  and  to  recommend 
those  most  suitable  for  aiding  the  planning  of  amphibious  operations.  This  was 
achieved  by  first  reading  relevant  literature  and  listing  the  different  types  of  models 
foxmd  within  it.  Then  their  relevant  attributes  (for  example:  the  equations  and 
assumptions  utilised  in  developing  them,  the  parameters  needed  to  utilise  them,  any 
vahdation  exercises  performed,  their  findings  and  their  accuracy)  were  analysed.  An 
important  observation  made  here  was  that  some  potentially  useful  numerical  models 
could  not  be  recommended  for  usage  to  the  RAN  as  they  required  parameters  which 
were  too  difficult  to  acquire. 
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Executive  Summary 

There  are  a  number  of  nearshore  characteristics  that  can  have  profound  effect  upon 
the  success  of  amphibious  operations.  Such  characteristics  include  wave 
amplitudes,  wave  propagation  angles,  longshore  and  cross-shore  velocities,  sand 
transport  rates  (an  indicator  of  turbidity)  and  beach  profile  characteristics  and  their 
temporal  changes.  Hence,  it  was  envisaged  that  a  set  of  reports  outlining  the  most 
appropriate  methods  (or  models)  for  calculating  such  parameters  would  be  of  great 
benefit  to  those  planning  amphibious  operations. 

In  July  1999,  General  Document  0214  was  produced  which  reviewed  nearshore 
wave  models  and  made  recommendations  as  to  which  were  suitable  to  be  used. 
The  aim  of  this  investigation  was  to  continue  this  work  and  produce  a  document 
which  outlined  the  different  types  of  nearshore  currents,  turbidity,  and  morphology 
models.  The  aim  was  to  then  review  their  attributes  (for  example:  (1)  the 
equations  and  assumptions  utilised,  (2)  parameters  needed  to  utilise  them,  (3)  any 
validation  exercises  performed  and  their  findings)  and  make  recommendations  to 
the  Australian  Defence  Force  (ADF)  as  to  which  models  were  most  suitable. 

An  important  observation  made  here  was  that  some  potentially  useful  numerical 
models  could  not  be  recommended  for  usage  to  the  RAN  as  they  required 
parameters  which  were  too  difficult  to  acquire. 
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1.  Introduction 


1.1  Stimulus  for  Report 

There  are  a  number  of  nearshore  characteristics  that  can  have  profound  effect  upon  the 
success  of  amphibious  operations,  or  military  operations  ocauring  within  the  coastal  zone. 
Such  characteristics  are  outlined  in  Del  Balzo,  Vodola  and  Beveridge  (1999)  and  from  these 
those  deemed  to  be  of  high  importance  included,  wave  amplitudes,  wave  propagation 
angles,  longshore  and  cross-shore  current  speeds,  sediment  transport  rates  (an  indicator  of 
turbidity)  and  morphology  characteristics  and  changes. 

To  aid  the  planning  of  amphibious  operations  it  was  envisaged  that  it  would  be  beneficial  if 
a  set  of  reports  was  written  outlining  the  most  appropriate  methods  (or  models)  for 
calculating  the  above  parameters. 

In  June  1999,  the  first  report  (namely  GD-0214  (Andrew,  1999))  was  produced  which 
outUned  the  most  current  nearshore  wave  models  and  made  recommendations  as  to  which 
were  most  applicable  to  the  needs  of  the  Australian  Defence  Force  (ADF).  This  report 
continues  on  from  where  GD-0214  finished  and  describes  those  models  which  forecast  those 
nearshore  characteristics  described  above  which  were  not  covered  in  GD-0124. 

This  DSTO  report,  summarising  the  most  appropriate  ways  of  predicting  or  calculating 
environmental  factors  pertinent  to  amphibious  operations,  is  also  important  to  otiier  DSTO 
tasks.  Most  noticeably,  the  models  cited  in  this  report  may  one  day  be  coded  into  the 
Amphibious  Warfare  Operational  Concept  System  Software  (AWOCSS)  being  developed  by 
Mr  Jamie  Watson  (DSTO  Task  #  98/165). 

1.2  Structure  of  Report 

In  this  report  pertinent  backgroimd  information  regarding  nearshore  processes  will  first  be 
outiined  (except  background  information  associated  with  surface  waves;  such  information 
was  outlined  in  GD-0214).  Then  the  different  types  of  nearshore  models  discovered  wiU  be 
summarised,  with  this  report  concluding  by  outlining  the  most  appropriate  models  for  the 
Australian  Defence  Force  (ADF)  based  upon  the  models  attributes  (assumptions  used, 
parameters  required,  validation  procedures  and  models  acctuacy). 

Note:  As  in  the  first  literature  review  (GD-0214),  it  will  be  asstimed  throughout  this  report  that 
surface  wind  conditions,  deep-water  wave  characteristics,  tidal  action,  shelf  current  regimes,  and 
static  bathymetry  characteristics  are  readily  available  for  any  nearshore  forecasting. 

1.3  Background  Information 

At  some  point  as  a  wave  propagates  towards  a  coastline,  the  drag  caused  b)^  the  wave 
"feeling"  the  bottom  wUl  cause  tiie  top  stages  of  the  wave  to  move  forward  at  a  greater  speed 
than  those  at  the  bottom  At  this  point  the  wave  wiU  become  hydrostaticaUy  tmstable  and 
break  with  such  a  position  marking  the  seaward  extent  of  the  surf  zone  and  being  strongly 
dependent  upon  the  depth. 
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As  at  any  location  there  is  normally  a  spectrum  of  wave  heights  (as  well  as  a  spectrum  of 
periods)  then  due  to  wave  breaking  being  strongly  dependent  upon  depth,  there  is  a  region 
in  which  all  the  waves  present  wiU  break.  This  region  is  part  of  the  surf  zone  and  for 
obvious  reasons  is  called  the  breaker  zone  (Fredsoe  and  Deigaard,  1992),  The  seaward  extent 
of  the  breaker  zone  is  where  the  waves  with  the  largest  height  break  with  the  smaller  waves 
breaking  further  inshore. 

When  waves  break,  their  momentum  is  transformed  into  "plunging",  vertical  motion,  as  well  as  a 
strong  horizontal  motion.  The  momentum  flux  of  this  horizontal  motion  has  been  given  the 
term  radiation  stress  (Longuet-Higgins  and  Stewart,  1964)  with  this  body  of  water  having  the 
ability  to  continue  aU  the  way  up  until  it  has  reached  the  top  stages  of  the  beach  face.  Here,  it  is 
termed  swash  with  the  term  backwash  applied  to  the  swash  that  runs  back  down  the  beach  face. 

In  the  surf  zone  there  is  a  group  of  waves  which  are  of  a  very  low  frequency  (0.005Hz  - 
0.05Hz)  and  whose  presence  can  influence  the  formation  of  longshore  currents,  rips  and 
longshore  bars.  These  waves  are  collectively  called  infragravity  waves  because  of  their 
location  in  the  power  spectra  of  surface  gravity  waves.  However,  this  collective  name  is 
somewhat  of  a  misnomer  as  they  are  still  gravity  waves  (Van  Rijn,  1998). 

In  the  ocean  there  are  three  types  of  infragravity  waves:  Boxmd  Long  Waves,  Leaky  Waves 
and  Edge  Waves  (Short,  1999).  The  last  two  of  these  are  important  for  the  formation  of  rip 
currents  and  multiple  sand  bars. 

Boimd  Long  Waves  are  associated  with  the  wave  packets  or  wave  groups.  They  are  a  lowering 
and  increase  in  the  mean  water  level  and  arise  from  the  fact  tiiat  the  large  waves  in  a  wave  group 
have  a  higher  radiation  stress  than  the  smaller  waves  in  the  wave  groups.  Hence,  horizontal 
radiation  stress  gradients  are  setup  in  which  it  has  been  theoretically  demonstrated  by  Longuet- 
Higgins  and  Stewart  (1964)  that  fluid  is  expelled  from  ti\e  high  waves  in  the  wave  groups  to  the 
lower  waves  in  the  wave  groups.  This  causes  a  reduction  in  the  mean  water  level  at  the  position 
of  the  large  waves  and  an  increase  in  the  mean  water  level  at  the  position  of  the  lower  waves  in 
the  wave  group.  Hence  the  boxmd  long  wave  is  180  degrees  out  of  phase  with  the  wave  groups 
(Short,  1999).  Figure#!  shows  an  example  of  a  Bound  Long  Wave. 


FigureUl  A  diagram  showing  the  position  of  a  bound  long  xuave  with  respect  to  the 
wave  groups.  Tahenfrom  Short  (1999). 
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Leaky  Waves  are  infragravity  waves  which  have  propagated  into  the  surf  zone  from  deeper 
water  but  have  then  been  reflected  back  off  the  shoreline  and  back  out  to  sea.  Hence  energy  is 
"Leaking"  out  of  the  surf  zone.  Their  wave  height  is  largest  at  the  shoreline,  is  constant  in  the 
longshore  direction  and  varies  in  the  offshore  direction  by  a  succession  of  alternating  maxima 
and  minima  of  diminishing  magnitude  (see  Figure#2).  Leaky  waves  form  when  Bound  Long 
Waves  approach  a  beach  nearly  normally  and  are  released  from  the  wave  groups  when  the  short 
period,  incident  surface  waves  break  and  the  wave  groupness  is  destroyed  (List,  1992).  The 
released  bound  long  wave  then  propagates  shoreward  as  a  free  wave.  At  the  shoreline  it  is 
reflected  back  out  to  sea  with  this  wave  then  having  the  ability  to  interact  with  those  propagating 
into  the  surf  zone  and  thus  generate  standing  leaky  waves  (Short,  1999). 


Figure#!  The  surface  elevation  signature  of  leaky  infragravity  waves  and  edge  waves  in  the 
surf  zone.  Taken  from  Short  (1999). 

Edge  waves  are  very  low  frequency  waves  that  travel  in  the  alongshore  direction  being 
trapped  within  the  stuf  zone  by  depth  induced  refraction.  They  have  their  largest 
amplitudes  at  the  shoreline,  with  the  amplitude  decaying  exponentially  in  the  offshore 
direction  and  sinusoidally  in  the  longshore  Section.  Figure#2  shows  an  example  of  a  simple 
edge  wave.  Edge  waves  of  the  same  frequency,  amplitude,  and  mode  but  travelling  in 
opposite  directions  can  constructively  interfere  to  produce  standing  edge  waves.  Such  edge 
waves  are  important  to  the  formation  of  multiple  sand  bars  and  have  the  same  frequency  as 
the  surface  incident  gravity  waves  (Short,  1999). 

There  are  two  formation  mechanisms  for  edge  waves.  For  steep,  reflective  beaches 
subharmonic  edge  waves  (those  having  periods  twice  that  of  the  incident  waves)  of  zero 
mode  are  produced  by  first  incident  waves  being  reflected  off  the  shoreline  to  produce 
standing  waves  which,  because  of  a  non-linear  transfer,  are  transformed  into  edge  waves. 
For  shallow,  dissipative  beaches,  Gallagher  (1971)  demonstrated  that  edge  waves  could  form 
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as  a  result  of  triad  non-linear  wave  to  wave  interaction  in  which  two  obliquely  incident 
bound  long  waves  propagating  towards  the  shoreline  could  superimpose  and  generate  a 
third  wave  (the  edge  wave)  of  the  same  frequency  as  the  incident  waves  (i.e.  synchronous 
edge  waves)  with  the  edge  waves  not  being  isolated  to  zeroth  mode  (Short,  1999). 

The  process  of  waves  approaching  a  coastline,  shoaling,  breaking,  and  running  up  the  beach 
face  is  relentless  and  has  the  ability  to  generate  longshore  velocities  between  the  breaker 
zone  and  the  swash  zone  which  eventually  feed  rip  currents. 

To  understand  how  longshore  currents  and  rips  form  in  the  ocean,  it  first  must  be  vmderstood 
that  the  decrease  in  the  onshore  radiation  stress  due  to  bottom  friction  is  directly  balanced  by  an 
increase  in  the  wave  setup.  Figure#3  shows  this  in  detail.  It  can  be  deduced  from  this  figtire  that 
if  there  is  a  longshore  variation  in  the  wave  setup  (caused  by  a  longshore  variation  in  the  wave 
height)  then  a  longshore  current  is  produced  in  which  water  will  nm  down  the  pressure  gradient 
imtil  it  converges  with  other  longshore  currents  generated  the  same  way  with  the  longshore 
currents  then  turning  seaward  and  producing  rips.  Figure#4  shows  an  example  of  this  scenario. 


Figure#  3  An  example  of  how  two  waves  with  different  offshore  xvave  heights  produce  different 
degrees  of  wave  setup  along  a  shoreline.  (Taken  from  Komar  (1998)). 

From  this  analysis  it  should  be  deduced  that  the  formation  of  rip  cells  is  solely  related  to 
longshore  variations  in  wave  heights.  In  the  ocean  wave  heights  can  vary  in  the  longshore 
direction  by  two  independent  mechanisms  with  each  mechanism  producing  markedly 
different  rip  cxirrents.  One  mechanism  is  whereby  the  beach  morphology  is  non-'variant  in 
the  longshore  direction  and  has  edge  waves  present  which  are  either  standing  or 
progressive.  As  the  edge  waves  have  small  wave  heights  which  vary  sinusoidally  along  the 
beach  then  this  causes  the  surface  wave  heights  and  wave  setup  to  vary  sinusoidally  along 
the  beach.  This  intum  produces  an  alongshore  pressure  gradient  which  produces  longshore 
currents  which  flows  seaward  as  rips  at  tiie  nodes  of  the  edge  waves  (Komar,  1998). 
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FigureM  A  schematic  picture  showing  hoio  xoaves  of  differing  heights  along  a  shoreline 
generate  rip  currents.  Taken  from  Komar  (1998). 

The  second  mechanism  is  where  rips  form  by  wave  refraction  resulting  from  bathymertiy 
variations  (headlands  or  offshore  canyons  or  ridges)  or  wave  diffraction  due  to  offshore 
islands.  The  wave  refraction  or  diffraction  causes  wave  convergence  and  divergence  zones 
along  the  coastline  which  results  in  the  wave  heights  to  vary  in  the  longshore  direction.  Like 
the  edge  wave  formation,  this  longshore  variation  in  the  wave  height  of  the  incident  waves  is 
translated  into  a  longshore  variation  in  the  wave  setup  which  leads  to  rip  formation. 
Figme#5  shows  the  location  of  two  permanent  rips  off  CaUfomia  which  originate  due  to  the 
presence  off  two  offshore  canyons.  For  the  rips  forming  as  a  result  of  headlands  these  are 
called  topographic  rips  and  can  sometimes  be  much  larger  in  volume  than  those  caused  by 
the  presence  of  edge  waves.  Such  large  topographic  rips  occur  under  storm  conditions  and 
are  called  mega-rips  (Short,  1999). 


Figure#5  The  location  of  two  permanent  rips  near  Scripps  Institution  of  Oceanography  (SIO), 
California,  and  their  close  proximity  to  two  offslwre  canyons.  Wave  heights  along  the 
shoreline  are  as  shown.  Talwnfrom  Komar  (1998). 
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As  far  as  the  beach  profile  is  concerned,  erosion  first  occurs  with  beach  material  being  placed 
into  suspension  by  the  strong  plunging  vertical  motion  accompanying  breaking  waves.  This 
plunging  motion  creates  an  offshore  transport  over  the  sea  bed  which  when  it  converges 
with  the  slight  onshore  transport  associated  with  non-hnear  (unbroken)  surface  waves 
results  in  the  generation  of  a  longshore  sand  bar  (Short,  1999).  Figure#6  is  a  schematic 
diagram  showing  such  a  formation  process.  Erosion  of  beach  material  can  then  continue  as  a 
result  of  the  "wall"  of  water  formed  by  the  breaking  waves  propagating  towards  the 
shoreline  and  either  moving  material  along  the  bed  (termed  bed  load  transport)  and/or  by 
placing  more  material  into  suspension  (termed  suspended  load  transport).  Figure#7  is  a 
schematic  diagram  illustrating  the  concept  of  bed  load  and  suspended  sediment  load. 


Figure#6  A  schematic  diagram  showing  how  a  longshore  bar  forms  from  the  convergence  of  sand 
transported  offshore  hy  breaking  ivaves  and  the  onslwre  sand  transported  associated  with 
onshore  wave  movement,  qsx  is  the  offshore  suspended  sediment  transport  rate  xvith  arrows 
lengths  indicating  relative  strengths.  (Modified  from  Fredsoe  and  Deigaard,  1992)). 


Figure#?  A  schematic  diagram  illustrating  tlxe  difference  between  bed  transport  load  and  suspended 
transport  load.  (Taken  from  Soulsby,  1997). 
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Bed  load  and  suspended  load  can  occur  anywhere  from  the  breaker  zone  to  the  swash  zone, 
with  their  magnitudes  depending  upon  the  velocity  of  the  broken  waves  and  the  particle  size 
distributions  of  the  beach  material.  Figure#8  shows  the  dependence  of  erosion  and  accretion 
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upon  these  two  parameters.  Bed  load  magnitudes  adjust  relatively  quickly  to  changing  wave 
and  current  conditions:  Suspended  loads,  however,  take  longer  to  respond.  This  is  because  it 
takes  a  finite  time  and  depth  for  new  sediment  on  the  bed  to  be  diffused  up  through  the  flow 
and  for  sediment  in  suspension  to  settle  (Soulsby,  1997).  Once  sand  and  beach  profile 
material  has  been  placed  into  suspension,  it  can  then  be  transported  to  other  locations  due  to 
the  action  of  any  longshore  currents. 


Figure^S  The  dependence  of  erosion  and  deposition  of  weathered  beach  sediment  upon  liorizontal 
velocity  and  partick  sizes.  (Taken  from  Bearman,  1989). 

In  addition  to  single  longshore  sand  bars  forming  as  a  result  of  the  plunging  or  vertical  nature 
of  broken  waves,  multiple  sand  bars  can  form  as  a  result  of  standing  infragravity  waves  (edge 
or  leaky)  waves.  In  such  a  situation  the  standing  edge  waves  form  as  a  result  of  the  interaction 
of  two  progressive  edge  waves  (which  may  be  propagating  in  the  same  or  in  opposite 
directions)  with  standing  leaky  waves  caused  by  the  interaction  of  leaky  waves  and  incoming 
bound  long  waves.  Such  standing  infragravity  waves  induce  horizontally  segregated  drift 
velocities  close  to  the  bed,  with  the  drift  converging  at  the  nodes  at  the  bottom  of  the  boimdary 
layers,  while  the  convergence  is  towards  the  antinodes  somewhat  higher  in  this  layer.  Figure#9 
shows  these  drift  velocity  peculiarities  for  standing  infragravity  waves.  Depending  upon 
whether  the  mode  of  sediment  transport  is  by  bed  load  or  suspended  load  will  determine 
where  any  longshore  bars  will  form.  However,  as  longshore  bars  form  in  high  energy  surf 
zones,  where  suspended  load  dominates,  then  longshore  bars  will  more  often  form  at  tiie 
antinodal  positions  of  tiie  standing  waves  (Short,  1999). 
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Figure#9  The  horizontal  velocity  (u)  orbital  characteristics  beneath  a  standing  infragravity  wave 
and  the  location  of  the  resulting  sand  bars  wlwn  the  sediment  movement  is 
predominantly  in  suspended  form  or  in  bedload  form.  Taken  from  Short  (1999). 
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2.  Longshore  Velocity  Predictions 

2.1  Monotonically  Sloping  Beaches 

The  very  first  models  of  longshore  velocity,  catised  by  obliquely  breaking  waves  on 
monotonically  sloping  beaches  which  had  no  rips,  were  proposed  by  Putman,  Munk  and  Traylor 
(1949),  Inman  and  Quinn  (1951)  and  Eagleson  (1965)  (Horikawa,  1988).  These  models  have  been 
reviewed  by  Galvin  (1967)  and  Komar  (1976),  and  have  been  largely  abandoned  (Komar  and 
Oltman-Shay,  1990).  Figure#10  shows  the  circulation  scheme  which  exists  under  such  conditions 
and  shows  that  the  largest  longshore  velocities  occur  close  to  the  mid-surfzone  position. 


FigureUlO  A  schematic  diagram  showing  the  horizontal  circulation  which  forms  when  waves 
break  obliquely  to  the  shoreline.  Adapted  from  Komar  (1998). 

The  first  meaningful  longshore  velocity  models  were  simultaneously  developed  by  Longuet- 
Higgins  (1970a,b),  Bowen  (1969)  and  Thornton  (1970).  Such  models  assumed  the  waves 
approaching  any  beach  were  monochromatic  and  stationary  in  time,  that  the  beach  profile 
had  a  uniform  slope,  and  that  the  bottom  contours  were  straight  and  parallel.  To  aid 
obtaining  a  solution  for  the  longshore  velocity  these  models  also  assumed  that  the  waves  are 
breaking  roughly  parallel  to  the  shoreline  (cos  ttb  =  1 )  (where  is  the  wave  approach 


angle  at  breaking)  and  that  the  maximum  velocity  of  the  wave  orbitals  (Uj^)  are  large 
(  <  1 )  relative  to  the  longshore  velocity  being  predicted  (v).  This  last  assumption  is 

only  valid  when  (tti,  <  24°).  The  major  difference  in  these  models  was  the  way  in  which  the 
horizontal  eddy  viscosity  coefficient  was  formulated  (Komar  and  Oltman-Shay,  1990). 


The  model  by  Longuet-Higgins  (1970a,b)  (equation#!)  was  used  extensively  and  was 
recommended  in  The  Shore  Protection  Manual  (SPM)  (CERC,  1984).  In  supporting  their 
recommendation  in  the  SPM,  it  was  quoted  that  the  model  fitted  the  data  by  Putnam,  Munk 
and  Traylor  (1949)  and  the  laboratory  data  by  Galvin  and  Eagleson  (1965).  Komar  (1998) 
reviewing  this  assertion  disagreed  and  reverberated  such  views  using  more  data  sets. 
Figure#!!  shows  the  comparisons  between  the  equation  derived  by  Longuet-Higgins 
(1970a,b),  the  data  sets  by  Putnam,  Munk  and  Traylor  (1949)  and  the  Galvin  and  Eagleson 
(1965),  and  other  data  sets  Komar  (1998)  utilised.  It  is  strongly  believed  that  the  discrepancy 
between  the  model  by  Longuet-Fliggins  (1970a,b)  and  insitu  and  laboratory  data  is  the 
incorporation  of  the  beach  slope. 
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V  =  0.207^gHb  sin(2ab) 


(1) 


where: 

1.  V  is  the  longshore  velocity  (ms-^)  at  the  midsurf  position  (half  way  between  the  shorehne 
and  breakpoint) 

2.  Hb  is  the  wave  height  at  breaking. 

3.  S  is  the  beach  slope. 

4.  g  is  the  acceleration  due  to  gravity  =9.8ms-2. 


FigureUll  A  scatter  plot  showing  the  model  of  Longuet-Higgins(1970a,b)  and  the  data  sets  used 
to  test  it.  Taken  from  Komar  (1998). 

For  beaches  of  constant  slope  under  the  action  of  monochromatic  oblique  waves  the 
longshore  velocity  formula  developed  by  Komar  and  Inman  (1970)  for  the  mid  surf  position 
provides  the  best  predictions  (Komar  and  Oltman-Shay,  1990).  In  testing  the  model  by 
Komar  and  Inman  (1970)  on  beaches  which  there  were  no  rips  and  only  longshore  velocities, 
it  was  observed  that  equation#2  was  applicable  up  to  large  breaking  angles  of  45  degrees 
and  not  the  ceiling  of  24  degrees  as  inferred  by  one  of  the  assumptions.  The  mid  surf 
position  is  the  location  where  the  largest  longshore  velocities  normally  occur  with  the 
equation  developed  by  Komar  and  Inman  (1970)  as  follows: 


V  =  0.585^gHb  sin(2aj 


(2) 


where: 

1 .  all  the  parameters  have  been  previously  explained 

Recently,  a  model  has  been  formulated  which  takes  into  accoimt  that  waves  are  irregular 
with  well  defined  probability  density  functions  (but  which  are  still  assumed  to  be  narrow 
banded  in  frequency  and  direction).  The  most  successful  of  these  (Holman,  1995)  was 
formulated  by  Thornton  and  Guza  (1986).  In  this  formulation  a  model  is  developed  by 
balancing  the  on-offshore  gradient  of  the  wave  induced  longshore  momentum  flux  (or 
radiation  stress)  against  the  longshore  bottom  stress  and  Reynolds  stresses.  This  model 
assumes  that; 

•  the  waves  conditions  are  stationary. 
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•  that  the  bottom  contours  are  monotonic  in  the  offshore  direction  (although  the  model  can 
be  extended  to  non-monotonic  beaches  by  solving  the  equations  iteratively),  and  are 
straight  and  parallel  to  each  other. 

•  Incident  wave  angle  is  small, 

•  And  like  the  monochromatic  wave  model  of  Komar  and  Inman  (1970)  assumes  that  the 
maximum  velocity  of  the  wave  orbitals  (u^)  are  large  (V  <1)  relative  to  the 

longshore  velocity  being  predicted  (v). 

In  applying  the  above  assumptions  together  with  a  linear  bottom  shear  stress  and  Snells  law, 
Thornton  and  Guza  (1983)  derived  equation#3.  This  equation  is  the  general  solution  for  beaches 
of  an  arbitrary  offshore  beach  profile,  although  before  it  can  be  used  equation#4  needs  to  be 
solved  for  Hims  by  numerical  integration.  Thornton  and  Guza  (1983)  then  simplified  equation#3 
by  considering  only  beaches  which  were  monotonic  in  tire  offshore  direction.  In  doing  so  they 
derived  an  aiialytical  solution  for  the  longshore  velocity  on-offshore  profile  (equation#5)  with 
both  the  analytical  model  and  general  beach  profile  solution  only  applicable  to  relatively  shallow 
waters  (where  the  ratio  of  depth  to  wavelength  is  <  1/20). 


3  B^fpVg  sina„ 

4  c,/  C„ 


(3) 


where: 

1.  Subscript  zero  refers  to  offshore  conditions  (but  the  shallow  water  assumption  should 
still  be  valid). 

2.  fp  is  the  peak  wave  frequency. 

3.  Y  is  a  constant  relating  the  water  depth  to  breaking  wave  height. 

4.  Co  is  the  wave  celerity  or  phase  speed  in  deep  water 

5.  h  is  the  water  depth. 

6.  Cf  is  the  friction  coefficient. 


i?P8H-C,cosa  =  ^V5^g;^f,HL 


(4) 
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where: 


In 

Y"  tanP 

B^fp 

2.  B  is  a  breaker  index  parameter  and  accoimts  for  the  different  breaker  types.  This 
parameter  can  be  visualised  as  being  a  function  of  the  amount  foam  on  the  face  of  each  bore 
(Thornton  and  Guza,  1983)  and  is  expected  to  be  less  than  unity  for  spilling  breakers  and 
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near  unity  for  plunging  breakers.  Unfortunately  this  parameters  needs  to  be  calculated 
using  insitu  wave  height  data. 

The  drawback  of  this  analytical  model  (equation#5)  is  that  too  many  parameters  need  to  be 
calculated  from  in-situ  data  in  order  for  this  model  to  be  applied.  Such  parameters  are  the 

•  breaker  index  B, 

•  the  parameter,  Y , 

•  offshore  wave  characteristics  of  incident  wave  height,  frequency  and  propagation 
direction, 

•  surf  zone  bathymetry  in  the  offshore  direction, 

•  coastline  orientation  needs  to  be  calculated  using  hydrographic  data. 

The  other  drawback  from  an  operational  perspective  is  that  the  longshore  velocity  actually 
needs  to  be  measured  at  a  number  of  points  so  that  the  longshore  velocity  cross-shore  profile 
can  be  predicted.  The  reason  for  this  is  that  the  bottom  friction  coefficient  Cf  cannot  be 
measured  in-situ  and  is  calculated  by  measuring  the  above  parameters  and  then  adjusting 
this  coefficient  tmtil  the  difference  between  longshore  velocity  predictions  and 
measurements  is  minimised. 

Thornton  and  Guza  (1986)  performed  such  a  task  using  longshore  velocity  and  wave  height 
data  collected  at  Leadbetter  Beach,  Santa  Barbara,  California  from  the  30*  of  January  to  the 
23’^^  of  February  1980.  This  beach  was  chosen  as  the  wave  heights  are  nearly  homogeneous 
in  the  longshore  direction  and  the  bottom  contours  are  straight  and  parallel.  Although 
longshore  velocity  data  sets  did  exist  prior  to  this  publication  for  similar  beaches  as 
Leadbetter  Beach,  their  on-offshore  instrument  spacing  was  not  dense  enough  (compared  to 
the  wavelength)  to  be  useful  with  the  main  reason  for  obtaining  this  data  set  was  to 
adequately  measure  the  wave  directional  properties.  Figure#12  shows  the  position  of  the 
instruments  in  water  depth  that  were  greater  than  3m. 


Figurettll  A  schematic  diagram  shoiving  the  location  ofxoave  and  longshore  velocity  sensors  deployed 
at  Leadbetter  beach,  California  betzoeen  30‘>^  January  and  23''^  February  1980.  The  mean 
wave  period  was  14seconds,  luith  the  surface  elevation  of  such  a  wave  shown  to  scale. 
Taken  from  Thornton  and  Guza  (1986). 
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For  Leadbetter  beach,  Thornton  and  Guza  (1986)  analysed  the  data  available  from  the  30*  of 
January  to  the  23"^^  of  February  1980  and  first  restricted  themselves  to  data  from  the  2"^  to  6* 
of  February.  Such  a  restriction  was  because  this  time  interval  was  when  the  model 
assumption  of  narrow  banded  waves  in  frequency  and  direction  was  valid.  During  this 
restricted  time  interval  the  model  assumption  of  planar  beach  offshore  profiles  was  nearly 
valid  and  most  valid  on  February  4*  (see  Figure#13). 


Figure#13  On-offsJwre  beach  profiles  and  mean  profiles  of  Leadbetter  beach  from  the  2"“^  to  6*  oj 
February  1980.  Note  how  they  are  all  nearly  planar  (small  single  bars)  with  the  most 
planar  profile  occurring  on  the  4*  of  February.  Taken  from  Thornton  and  Guza  (1986). 

Using  offshore  and  inshore  wave  characteristics  on  the  4*  of  February,  the  wave  parameters 
B  and  the  parameter  Y  were  calculated.  Then  the  beach  orientation  was  deduced  using  the 
orientation  of  the  50m,  Om,  -50m  and  -100m  isobaths  referenced  to  the  Mean  Sea  Level. 
Lastly,  these  pieces  of  information  were  then  used  to  calculate  the  cross-shore  variation  of 
the  longshore  velocity  by  iteratively  solving  the  least  square  error  between  the  calculated 
and  measured  longshore  velocities  shoreward  of  the  breaker  line.  This  procedure  gave  the 
bottom  friction  coefficient  Cf  a  value  of  0.009.  The  calculated  cross-shore  profile  and  the 
longshore  velocity  data  used  to  generate  it  are  as  shown  in  Figure#14.  It  should  be  seen 
from  such  a  figure  lhat  the  model  appears  very  good,  however,  such  a  view  can  be 
misleading  as  to  achieve  such  a  comparison  the  bottom  friction  coefficient  was  optimally 
adjusted.  Although  a  reassuring  observation  is  that  the  value  reached  for  Cf  was  within  the 
range  of  what  has  been  observed  (0.005  to  0.3  (Shemdin  et-al.,  1978)). 
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Figureitl4  Predicted  and  measured  longshore  velocities  for  Leadbetter  Beach,  California  on  February 
4"’,  1980.  Taken  from  Thornton  and  Guza  (1986). 

2.2  Non-monotonically  Sloping  (or  barred)  Beaches. 

There  have  been  a  few  models  developed  which  attempt  to  estimate  the  strength  of  the 
longshore  current  on  a  barred  beach.  The  first  of  these  foimd  in  the  literature  was  by 
Deigaard,  Fredsoe  and  Hedegaard  (1986).  Although  such  a  model  appears  very  good  it  does 
have  the  common  problem  of  predicting  the  strongest  currents  directly  over  longshore  bars 
and  the  weakest  over  the  troughs.  A  characteristic  which  is  inconsistent  with  detailed 


FigurettlS  Calculated  longshore  velocities  (Vy)  for  a  barred  coastal  profile  for  tivo  offshore  wave 
characteristics.  Note  tlje  location  of  the  maximum  velocity  over  the  bar  crest  although  high 
density  data  has  it  located  in  troughs.  T  is  the  wave  period.  Taken  from  Fredsoe  and 
Deigaard  (1992). 

observations.  Figure#15  shows  this  characteristic. 

The  best  model  for  predicting  tbe  longshore  velocity  and  which  goes  someway  to  predicting 
the  strongest  longshore  currents  over  the  troughs  is  that  developed  by  Smith,  Larson  and 
Kraus  (1993)  (Holman,  1995).  This  model  was  developed  from  that  by  Larson  and  Kraus 
(1991)  by  including  a  Turbulent  Kinetic  Energy  term  in  the  radiation  stress  formulation.  Such 
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an  approach  shifts  the  location  of  the  maximiun  longshore  velocities  from  over  the  bars 
shoreward  towards  the  troughs.  The  model  developed  improved  the  agreement  between 
predictions  and  data  by  a  20-50%  reduction  in  least  squares  error  (Smith  et-al,  1993). 

2.3  Summary 

From  critically  reading  the  available  literature  it  was  concluded  that  the  best  model  for 
calculating  the  longshore  velocity  on  monotonicaUy  sloping  beaches  of  constant  slope  due  to 
the  action  of  monochromatic  oblique  waves  was  developed  by  Komar  and  Inman  (1970)  (as 
cited  by  Komar  and  Oltman-Shay,  1990).  If  the  irregular  nature  of  the  wave  heights  needs  to 
be  taken  into  account  then  the  model  by  Thornton  and  Guza  (1986)  will  yield  the  best  results 
with  the  model  by  Smith  et-al.  (1993)  (as  cited  by  Holman,  1995)  producing  the  best  residts 
for  non-monotonically  sloping  beaches. 

2.4  Recommendations 

By  analysing  the  formulations  of  the  above  models  and  considering  what  parameters  the 
ADF  is  likely  to  be  able  to  measure  or  predict,  the  only  model  which  can  be  recommended  is 
the  very  basic  and  least  accurate  model  by  Komar  and  Inman  (1970).  This  is  because  this  is 
the  only  model  for  which  aU  its  input  parameters  (for  example  surface  wave  characteristics) 
can  be  predicted:  The  model  by  Thornton  and  Giiza  (1986)  would  be  difficult  to  implement 
as  bottom  friction  needs  to  be  measured,  with  the  model  by  Smith  et-al.  (1993)  xmtenable  as 
turbulent  kinetic  energy  needs  to  be  measured. 
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3.  Two  Dimensional  Wave  Induced  Flow: 

3.1  Introduction 

It  can  be  easily  stated  that  no  numerical  models  exist  which  predict  the  strength  and 
direction  of  the  horizontal  current  in  the  surf  zone  (which  can  be  decomposed  into  longshore 
currents  and  any  rip  currents)  for  arbitrary  bathymetry  and  wave  conditions  as  a  function  of 
location  as  weU  as  time.  Such  a  model  would  need  as  input,  the  time  varying  surface  waves, 
as  well  as  high  resolution  bathymetry  input  of  the  surf  zone  and  the  location  of  any  edge 
waves  as  a  function  of  time.  The  position  of  edge  waves  as  functions  of  time  is  not  possible  to 
predict  with  simply  the  detection  of  edge  waves  very  difficult,  requiring  collecting  current 
data  in  the  energetic  surf  zone  from  numerous  sensors.  Hence  models  which  predict  the 
location  and  three  dimensional  structure  of  any  rips  do  not  exist  and  are  unlikely  to  for  at 
least  several  years. 

What  does  exist  however,  is  basic  quantitative  models  which  permit  the 

•  basic  characteristics  of  rips  associated  with  rip  channels  (formed  by  the  presence  of 
edge  waves)  to  be  quantified. 

•  longshore  spacing  of  edge  wave  induced  rips, 

•  formation  of  mega  rips  to  be  predicted, 

3.2  Simple  Models  of  Edge  Wave  Induced  Rips 

Rips  along  any  beach  can  be  sighted  from  the  air  or  from  the  beach  by  the  higher  turbidity 
associated  wiA  them  compared  to  neighbouring  waters.  Rips  are  also  likely  to  be  visually 
detected  in  the  surf  zone  by  the  decrease  in  the  wave  height  they  cause  for  the  incoming 
waves,  as  the  rip  extends  seaward  out  past  the  breaker  zone  (Short,  1999).  The  decrease  in 
the  wave  height  of  the  incoming  waves  because  of  a  the  rip  moving  in  opposite  direction  to 
the  waves  propagation,  is  due  to  the  principle  of  conservation  of  energy. 

The  horizontal  velocity  of  a  rip  as  it  extends  seaward  can  vary  greatly  from  very  low 
velocities  (0.1  to  0.5m/ s)  to  1  or  2  m/s.  The  strength  of  the  mean  rip  velocity  is  a  fimction  of 
incident  wave  height  as  well  as  the  phase  of  the  local  tides,  with  rip  velocities  increasing 
under  storm  conditions  and  also  increasing  vmder  low  tidal  conditions  due  to  the  larger 
wave  dissipation  over  the  more  exposed  longshore  bars.  Under  low  tide  levels  the  longshore 
currents  and  any  rips  they  feed  are  also  stronger  as  the  flow  is  constrained  to  flow  in  any 
formed  longshore  troughs  and  rip  channels.  The  strength  of  any  rips  as  they  extend  seaward 
has  also  been  observed  to  pulsate  with  periods  of  aroimd  30  seconds  to  several  minutes  with 
this  phenomena  attributed  to  far  infragravity  waves  (called  shear  waves)  (Short,  1999). 

The  longshore  spacing  between  non-topographic  rips  is  approximately  four  to  eight  times 
the  surf  zone  width  and  depends  on  the  beach  slope  and  wave  period;  the  lower  the  beach 
slope  the  wider  the  spacing  with  the  spacing  also  increasing  for  waves  of  longer  period.  On 
non-sheltered  beaches  rip  spacing  can  vary  from  100m  to  500m  although  commonly  rips  are 
150m  to  250m  apart  (Short,  1999). 
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If  rips  form  as  a  result  of  edge  waves,  and  hence  form  at  the  edge  wave  nodes,  then  the 
wavelength  of  the  rips  will  match  that  of  the  edge  waves.  If  the  distance  between  rips  needed 
to  be  forecasted  then  the  wavelength  of  the  edge  waves  could  be  used  as  a  proxy. 
Unfortunately  the  detection  of  edge  waves  is  difficult  and  time  cor^uming  (instrument 
deployment  and  data  analysis  is  approximately  several  days)  and  needs  several  longshore  or 
offshore  evenly  spaced  instruments  in  the  vigorous  surf  zone  (for  an  example  see  Oltman- 
Shay  (1998)).  ha  terms  of  predicting  the  edge  wavelength  (Le),  this  is  a  function  of  the  edge 
wave  period,  the  beach  slope  and  the  edge  wave  mode,  as  given  in  equation#6:  Again, 
parameters  which  need  to  be  measured.  Hence,  it  would  be  difficult  at  tiie  moment  for  the 
ADF  to  forecast  the  spacing  of  rip  currents.  If  the  spacing  between  rip  currents  is  important 
then  photographs  taken  from  low  flying  aircraft  would  yield  the  best  results. 


L.  =:^’g(2n+l)tanpl 


(«) 


where: 

1 .  Le  is  the  wavelength  of  the  synchronous  edge  wave 

2.  Te  is  the  edge  wave  period. 

3.  n  is  the  edge  wave  mode  number. 

4.  P  is  the  slope  of  the  surf  zone. 

3.3  Longshore  Bars  broLen  by  Regularly  Spaced  Rip  Channels 

3.3.1  Dalrymple  (1978) 

The  basic  model  by  Dalrymple  (1978)  was  developed  to  calculate  the  onshore  movement  of  water 
(q)  over  a  longshore  bar  and  subsequent  longshore  velocity  (V)  within  the  trough  due  to  normally 
incident  breal^g  waves.  The  bar  in  this  model  was  of  length  2L  and  was  dissected  by  evenly 
spaced  rip  channels.  Such  rips/rip  channels  form  due  to  Ae  presence  of  edge  waves  (either 
standing  in  the  longshore  direction  or  progressive).  The  driving  force  for  the  longshore  flow  in  the 
trough  and  subsequent  offshore  flow  in  the  rip  channels  is  longshore  variations  of  wave  setup, 
with  such  variations  being  induced  by  rip  channels’  effect  on  wave  refraction  and  because  wave 
breaking  in  the  deeper  rip  channels  is  less  intensive.  Figure#16  is  a  schematic  showing  the  surf 
zone  geometry  of  a  single  bar  of  length  2L  which  is  bordered  by  evenly  spaced  rip  channels. 


Figure#16  A  schematic  diagram  showing  the  surf  zone  geometry  of  the  model  developed  by 
Dalrxfmple  (1978).  See  text  ibr  explanation  of  variables. 
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The  model  developed  assumes: 

•  normally  incident  waves, 

•  the  width  of  rip  channels  are  large  enough  so  that  the  longshore  flow  shoreward  of  any 
two  neighbouring  bars  does  not  interfere. 

To  determine  the  magnitude  of  the  onshore  current  over  the  bar  Dairy mple  assumed  that  the 
loss  in  energy  was  equal  to  the  velocity  head  of  the  net  discharge  over  the  bar.  He  also 
assumed  that  the  bed  shear  was  not  important,  with  later  analysis  showing  that  such  an 
approach  was  justified.  The  magnitude  of  the  onshore  flow  (q)  is  given  by  equation#7a. 


q=D,^2d.m,-m) 


(7a) 


Where: 

1 .  q  is  the  onshore  flow  of  water, 

2.  Dc  is  the  water  depth  over  the  bar  crest, 

3.  ADq  is  the  setup  when  there  is  no  net  current  and  uniform  conditions 

4.  AD  is  the  actual  setup  when  there  is  an  onshore  current. 

To  calculate  the  longshore  velocity  Dahymple  used  the  longshore  projected  momentum 
equation  for  non-uniform  flow  as  well  as  the  continuity  equation.  An  assumption  utilised  in 
this  approach  was  that  the  wave  setup  is  small  compared  to  the  water  depth  in  the  trough, 
which  is  a  valid  assumption,  and  also  to  neglect  bed  shear  stress.  The  longshore  velocity 
equation  developed  is  as  shown. 


^  J 

sinH 


V2D,y 

A 


\ 

-tanl 


V2D,L 

A 


(7b) 


Where: 

1.  A  is  the  cross-sectional  area  of  the  trough. 

2.  is  simply  the  gravitational  constant  multipHed  by  the  wave  setup,  although  wiU 
be  denoted  by  the  velocity  Vo  as  it  has  the  same  units. 

To  illustrate  an  example  of  using  this  equation,  a  longshore  bar  of  200meters  in  length  was 
chosen,  which  had  a  depth  of  Imeter  at  its  crest  and  a  trough  of  area  lOOm^  (calculated  using 
the  geometry  shown  in  Figure#! 6).  An  offshore  wave  height  of  2  meters  with  a  steepness  of 
1/20  was  also  chosen.  Using  the  wave  steepness,  wave  height  and  depth  over  the  bar  crest 

of  Im,  Figure#!?  from  Dahymple  (1978)  was  used  to  deduce  the  magnitude  of  yJgAD^ 

(represented  as  Vo)  in  equationTb,  with  a  value  of  !.6m/s  reached.  Using  this  value 
equationTb  was  then  used  to  calculate  the  variation  of  the  longshore  velocity  along  the 
shoreline  with  Figure#!8  generated. 

From  Figure#!8  the  longshore  symmetry  expected  in  the  longshore  velocity  for  a  bar 
dissected  by  evenly  spaced  rip  channels  was  observed  with  the  magnitudes  predicted  in  the 
range  one  would  expect.  Although  no  test  could  be  found  for  this  model  the  results 
generated  in  Figure#!8  illustrate  that  such  a  model  appears  very  worthwhile. 
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Figure#17  Tlte  relationship  between  the  parameter  ^gAD^  (defined  as  Vo)  as  a  function  of  offshore 

wave  steepness  (Ho/Lo)  for  differing  depths  (Doc)  at  the  bar  crest.  Taken  from  Dalrymple 
(1978). 


Figure#18  Longshore  variation  in  tlie  longshore  velocity  in  the  trough  in  Figure#16  using 
Daln/mple's  model.  See  text  for  beach  characteristics. 

3.4  Simple  Model  of  Mega  Rips 

In  the  introduction  of  this  report  it  was  outlined  that  rips  can  form  due  to  the  presence  of 
headlands  (topographic  rips),  due  to  the  presence  of  edge  waves  or  due  to  unique  offshore  wave 
refraction  schemes.  For  the  first  formation  scheme,  a  simple  model  has  been  generated  which 
permits  the  impact  headlands  have  on  the  cell  circulation  caused  by  edge  wave  rips  to  be 
determined.  The  model,  which  does  not  take  into  accoxmt  tidal  effects,  was  developed -by 
Martens,  Williams  and  Cowell  (1999)  and  appropriately  uses  a  parameter  called  the  non- 
dimensional  emba5mient  scaling  parameter  (5  )  (NDESP)  (eqn.8).  This  parameter  can  be  used  to 
predict  if  topographic  rips  will  form  parallel  to  headlands  which  are  encapsulating  an 
embayment  but  it  should  be  noted  that  the  NDESP  cannot  be  used  to  infer  if  rips  will  form  by 
edge  waves  in  the  interior  of  embayments. 
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(8) 


where: 

1 .  Sa  is  the  shoreline  length  between  headlands  (see  Fig#19). 

2.  Cl  is  the  width  of  the  embayment  (see  Fig#19). 

The  non-dimensional  embayment  parameter  is  a  function  of  the  distance  between  two 
neighbouring  headlands  (along  the  shoreline  and  also  the  shortest  distance  through  the  surf 
zone)  and  also  the  wave  height  at  breaking.  Figure#19  is  a  schematic  diagram  showing  the 
parameters  used  to  calculate  the  non-dimensional  embayment  parameter.  It  has  been 
observed  that 


•  topographic  rips  occur  at  either  end  of  an  embayment  when  the  NDESP  is  less  than  8. 
Under  such  conditions  the  circulation  scheme  is  termed  cellular  with  the  coastal 
dimensions  and  wave  conditions  not  permitting  edge  wave  rips  to  form  in  the  interior  of 
the  embayment.  Under  storm  conditions  topographic  rips  have  stronger  velocities, 
extend  further  out  to  sea  and  are  hence  termed  mega  rips. 

•  embayments  have  topographic  rips  and  maybe  edge  wave  rips  when  the  NDESP  is 
between  8  and  19.  The  circulation  in  this  situation  is  called  transitional. 

•  an  embayment  will  have  no  topographic  rips  when  the  NDESP  is  greater  than  19.  It  may 
however  have  rips  forming  due  to  the  presence  of  evenly  spaced  progressive  or  standing 
edge  waves.  The  circulation  scheme  under  this  situation  is  called  the  normal  circulation. 


Figure#!  9  A  schematic  showing  the  definition  of  the  parameters  used  to  calculate  the  non- 
dimensional  scaling  parameter.  Modified  from  Short  (1999). 

In  addition  to  the  influence  the  wave  heights  and  embayment  characteristics  have  on  the 
formation  of  topographic  rips,  the  beachface  slope  also  has  an  affect.  This  is  because  the 
beach  slope  controls  the  amount  of  wave  energy  reflected  back  out  to  sea  and  hence  the 
ability  for  standing  edge  waves  to  form.  The  more  standing  edge  waves  have  the  ability  to 
form  the  more  a  normal  rip  cell  circulation  wiU  prevail  which  wiU  limit  the  growth  of 
topographic  rips. 
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The  slope  of  a  beachface  cannot  be  measured  once  and  used  for  all  wave  heights  and  wave 
periods  as  it  will  change  with  changing  wave  regimes  and  is  also  a  function  of  the  sediment 
making  up  the  beachface.  Hence  a  parameter  needed  to  be  developed  which  allowed  the 
beach  slope  to  qualitatively  be  predicted  knowing  the  sediment  type,  wave  height  and  wave 
period.  This  parameter  is  the  dimensionless  fall  velocity  (DFV)  (Q )  and  was  developed  by 
Wright  and  Short  (1984)  for  natural  beaches  by  extending  the  laboratory  work  of  Gomlay 
(1968).  The  DFV  is  a  function  of  the  breaking  wave  height  (Hb),  wave  period  (T)  and 
sediment  fall  velocity  (®s)  (see  Eqn#9)  and  it  was  observed  on  natural  beaches  that  when 
Q<1  beaches  were  reflective  (steep  (>3°)  and  with  no  longshore  bars),  when  Q>6  the 
beaches  were  dissipative  (flat  and  with  multiple  longshore  bars),  and  with  beaches  being 
intermediate  when  2<  <5  (with  one  or  two  longshore  bars). 


Due  to  the  influence  the  slope  of  the  beachface  has  on  the  rip  circulations.  Short  (1999) 
extended  the  work  of  Martens  et  al.  (1999)  to  illustrate  the  impact  it  has  on  the  cell 
circulations  within  the  surf  zone.  Figure#20  shows  the  influence  the  type  of  the  beach 
(predicted  using  Q)  has  upon  the  prevaiUng  cell  circulation  for  steep  beaches  and  for 
dissipative  beaches.  It  can  be  seen  from  this  figure  that  for  reflective  beaches  (Q  =1-5)  with 
headlands  a  long  way  apart  (5  >20)  the  rip  circtdations  are  solely  the  result  of  the  edge 
waves.  As  the  headlands  become  closer  together  (6  <8)  for  the  same  reflective  regime 
(Q=l-5)  the  headlands  are  too  close  together  to  allow  edge  waves  to  form.  Hence  under 
^ese  conditions  the  small  mass  flux  of  water  onto  the  beach  is  available  to  drive  topographic 
rips.  For  dissipative  beaches  ( Q  >6),  which  are  beaches  for  which  the  ability  of  edge  waves 
to  form  diminishes  greatly,  there  is  ample  mass  transport  onto  beaches  to  permit 
topographic  rips  to  form.  For  such  embayments  the  topographic  rips  may  be  present  parallel 
to  both  headlands  or  only  one  depending  on  the  wave  approach  angles. 

In  the  above  paragraphs  regarding  predicting  the  formation  of  topographic  rips  in 
embayments,  the  embaymentisation  parameter  is  a  parameter  which  could  be  relatively 
easily  calculated  by  the  ADF  using  topographic  charts  and  wave  heights  from  numerical 
models.  To  calculate  the  dimensionless  iaH  velocity  ( Q )  however,  requires  knowledge  of  the 
sediment  fall  velocity  s  which  would  be  difficult  to  measure  in  the  field  by  the  ADF  and  as 
such  needs  to  be  calculated. 

Initially  it  was  thought  that  to  calculate  the  settling  velocity  ( 0)s)  would,  Uke  measuring  it, 
be  difficult.  By  consulting  Soulsby  (1997)  it  was  observed  that  the  settling  velocity  is 
relatively  easy  to  calculate  with  the  best  method  for  calculating  it  being  the  formula 
developed  by  Soulsby  himself  (HaUermeier,  1981)  (see  equation#10a).  The  advantage  of 
equation#10a  is  that  the  only  parameter  which  needs  to  be  guessed  (as  it  would  be  difficult 
to  measure)  is  the  median  grain  size  (d).  Such  a  parameter  can  be  guessed  by  taken  a  sample 
of  the  sediment  and  visually  analysing  it  for  content.  The  other  parameters  of  the  settling 
velocity  which  can  be  taken  as  constant  are  the  grain  and  sea  water  density  ratios  (s)  and  the 
kinematic  viscosity  of  sea  water  ( V ): 
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Figure#20  A  schematic  showing  the  rip  circulation  schemes  ivhich  form  as  a  function  of  the 
dimensionless  embaymentisation  parameter  and  the  beach  type  (depictured  by  the 
dimensionless  fall  velocity  ( Q )).  Taken  from  Short  (1999). 

•  the  density  ratio  of  sea  water  to  sediment  can  be  taken  as  a  constant  as  it  can  be 
assumed  beaches  are  composed  of  sand  ( Ps  =  2650  kgm-3)  and  that  die  sea  water 

density  can  be  a  constant  1027kgm-3 .  Hence  (s)  in  equation#10b  can  be  taken  as  2.58 
(dimensionless) . 

•  The  kinematic  viscosity  changes  so  little  (see  Figure#21)  that  it  can  be  taken  as  10-^ 
mVi,  although  it  can  be  accurately  calculated  knowing  the  temperature  of  the  sea 
water  and  is  such  a  small  function  of  salinity  that  including  salinity  in  the 
calculations  is  deemed  not  necessary.  Figure#21  shows  the  dependence  of  the 
kinematic  viscosity  on  the  water  temperature  for  fresh  water  (zero  salinity)  and  for.  a 
salinity  of  35psu  (practical  salinity  units). 


CO,  =^y(l0.36"  +1.049D^)- 10.36 


J  (10a) 
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where: 

1 .  V  is  the  kinematic  viscosity  of  the  sea  water 

2.  d  is  the  median  diameter  of  the  sediment 

3.  D»  is  the  dimensionless  grain  size  and  given  by  equation#  10b 


3 


(10b) 


where; 
1.  s  is 


Figure#21  The  dependence  of  kinematic  viscosity  on  temperature  for  freshxoater  (zero  salinity)  and  a 
salinity  of  35ppt  (parts  per  thousand).  Taken  from  Soulsby  (1997). 


From  the  above  paragraphs,  it  was  concluded  that  the  model  by  Martens  et-al.  (1999), 
enabling  the  presence  of  topographic  rips  to  be  predicted,  is  a  relatively  simple  model  to 
apply  and  would  be  a  valuable  tool  for  the  ADF.  The  two  components  of  this  model,  namely 
the  embaymentisation  parameter  and  the  dimensionless  fall  velocity,  are  also  relatively 
simple  to  calculate  with  all  that  is  needed  being  topographic  charts  and  wave  heights  for  the 
embaymentisation  parameter,  and  sediment  characteristics  for  the  dimensionless  fall 
velocity. 

3.5  Summary 

From  reading  die  available  literature  it  is  concluded  that  there  are  no  numerical  forecasting 
models  available  which  can  predict  the  strength  of  the  horizontal  current  in  the  surf  zone 
due  to  gradients  in  the  radiation  stress  as  a  function  of  time  and  location  and  that  such  a 
model  is  unlikely  to  be  available  for  a  long  time  (more  than  5  years).  What  has  been  shown 
are  the  magnitudes  of  basic  rip  characteristics  (velocities  and  longshore  spacings),  as  weU  as 
models  which; 

•  predict  the  longshore  spacing  of  rips  caused  by  the  presence  of  edge  waves, 

•  predict  the  rip  characteristics  for  longshore  bars  which  are  evenly  dissected  by  rip 
channels 

•  predict  the  ability  of  topographic  rips  to  form  in  embayments. 
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3.6  Recommendations 

From  reading  the  relevant  scientific  literature  it  can  be  stated  that  the  broad  scale 
characteristics  of  rips  (their  strengths  and  longshore  spacings)  are  known  and  that  such 
values  can  be  used  by  the  ADF.  The  scientific  literature  also  yielded  numerical  models  of  rip 
characteristics,  as  well  as  their  feeder  currents.  By  comparing  the  input  required  by  these 
models  with  the  capabilities  of  the  ADF  it  was  concluded  that  the  models  which  can  be 
recommended  to  the  ADF  include: 

•  The  model  by  Dalrymple  (1978),  which  permits  the  longshore  velocity  longshore 
variation  behind  a  single  bar,  which  is  evenly  dissected  by  rip  channels. 

•  The  model  by  Martens  et-al.  (1999)  which  permits  the  formation  of  topographic  rips 
in  embayments  to  be  predicted. 
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4.  Longshore  Sediment  Transport: 

4.1  Introduction 

For  military  operations  knowing  the  visibility  of  the  water  column  at  any  depth  is  a  very 
worthwhile  parameter  to  be  able  to  forecast.  Such  a  parameter,  for  example,  is  helpful  to 
clearance  divers  conducting  mine  warfare  operations.  One  parameter  of  visibility  is  the 
amoxmt  of  sediment  within  the  water  column.  In  Section#5  of  this  General  Document, 
models  which  permit  the  sediment  concentrations  as  a  function  of  depth  to  be  predicted  will 
be  outlined.  However,  before  this  takes  place  it  was  considered  that  as  a  precursor,  models 
which  permit  the  total  longshore  sediment  transport  to  be  predicted  within  the  surf  zone  due 
to  the  action  of  breaking  waves  (or  longshore  radiations  stress  gradients)  should  be  outlined. 
Such  models  are  not  directly  applicable  to  any  facet  of  military  operations  and  have  been 
summarised  here  only  for  completeness. 

4.2  Longshore  Sediment  Transport  due  to  Broken  Waves 

One  of  the  most  widely  used  methods  for  calculating  the  total  longshore  sand  transport 
(non-cohesive  sediment)  across  the  surf  zone  was  developed  by  CERC  (Coastal  Engineering 
Research  Centre)  (Soulsby,  1997).  This  formtda,  however,  can  only  be  appUed  to  beaches 
with  the  following  characteristics: 

1.  Fine  to  medium  sediment  grain  sizes  (diameter  <  0.5mm) 

2.  Beaches  with  no  tidal  currents  (thus  inhibiting  its  applicability  to  the  northern  Austrahan 
coastal  zone). 

3.  Low  beach  slope  <  1:100 

4.  Straight  beaches 

The  CERC  formula  has  been  superseded  by  the  Queens  formxda 
(http://www.ihe.nl/he/dicea/sed41/sed4102.htm).  This  is  because  Schoones  (1996) 
comparisons  of  tiie  Queens  formula  and  others  against  in-situ  data  revealed  that  it  is  the  best 
cxurently  available  (Holman,  1995).  Schoonees  also  highlighted  the  fact  that  the  Queens 
formula  is  a  lot  more  applicable  than  the  CERC  formula. 

The  Queens  formula  was  derived  from  dimensional  analysis  and  is  expressed  as  follows. 


(11) 


Where: 

(1)  Q  is  the  longshore  semd  transport  rate  in  Kg/  s 

(2)  K  is  a  dimensionless  constant  which  can  be  determined  from  field  data,  although 
Kamphuis  (1986, 1990),  using  an  extensive  data  set,  observed  a  value  of  1.3*10-3. 

(3)  is  the  sea  water  density 

(4)  H  is  the  mean  water  depth 

(5)  T  is  the  peak  period  of  the  spectrum  of  waves 
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(6)  dso  is  the  median  grain  size 

Despite  the  success  the  Queens  formula  has  had,  it  still  has  inhereted  restrictions.  Similar  to 
the  CERC  formula,  the  Queens  formula  should  only  be  applied  to  beaches  with  the  following 
characteristics: 

1.  no  tidal  currents  (which,  like  the  CERC  formula,  inhibits  its  applicability  to  the  northern 
Australian  coasthne). 

2.  straight  beaches  with  no  groynes  or  breakwaters. 

3.  plane  (monotonicaUy)  sloping  beaches  with  no  lonsghore  bars. 

4.3  Longshore  Sediment  Tranport  due  to  Waves  and  Currents 

As  specified  above,  the  CERC  and  Queens  formulas  should  only  be  applied  to  beaches  where 
the  horizontal  speed  of  any  external  currents  (ie  Ekman,  tidal,  geostrophic)  is  very  low 
(typically  <  20%  of  the  longshore  velocity  induced  by  waves)  Soulsby  (1997).  Hence,  for  the 
strong  tidal  regimes  of  the  northern  Australian  coastline  (see  Fig#22)  other  formulations 
needed  to  be  found. 


Figure#22  The  tidal  regimes  encountered  around  the  Australian  coastline.  Note  how  the  northern 
Australian  coastline  has  lar^e  tidal  fluctuations.  Modified  from  Harris  (1994). 

From  Soulsby  (1997),  four  models  are  described  for  nearshore  environments  which  may  have 
strong  external  currents.  The  first  was  by  Badard  (1981)  and  was  developed  using  the 
"energetics"  arguments  of  Bagnold,  Inman  and  Bowen.  Such  an  argument  can  be 
summarised  as  saying  that  the  work  done  in  transporting  sediment  is  a  fixed  proportion  of 
the  total  energy  dissipated  by  the  flow. 

The  model  formulated  by  Badard  was  for  bed  and  suspended  load  and  can  include  the 
effects  of  wave  asymmetry  and  bed  slopes.  Applying  the  model  however,  to  surf  zones  with 
rippled  seabeds  can  result  in  unrealistic  sediment  transport  predictions  with  the  model  also 
only  being  applicable  when  low  external  currents  exist. 
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For  bed  and  suspended  load  over  rippled  beds,  the  models  by  Soulsby-Van  Rijn  (see  Soulsby 
(1997)  page  183)  or  Van  Rijn  (1989)  should  be  used.  The  model  by  Sotilsby  was  developed  by 
applying  the  principles  of  Grass  (1981)  to  the  current  induced  transport  model  by  Van  Rijn 
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(1984)  and  modifying  it  to  include  a  sediment  threshold  term  and  a  bed  slope  term  (Soulsby, 
1997).  The  model  Van  Rijn  (1989)  was  developed  by  adapting  his  current  alone  formula 
(namely  Van  Rijn  (1984))  to  include  an  analytical,  semi-empirical,  model  of  sediment 
diffusion  through  the  wave  boundary  layer.  It  is  quite  a  complicated  model  and  calculates 
both  longshore  velocities  as  weU  as  suspended  loads  and  bed  loads.  Comparisons  of  the 
model  results  with  in-situ  data  are  quite  good  (Soulsby,  1997). 

The  fourth  model  is  the  Sediment  Transport  at  a  Point  model  (STP)  and  hke  the  Van  Rijn 
model  is  an  extensive  computer  based  model  which  calculates  longshore  velocities  as  well  as 
longshore  sediment  transports.  It  can  be  purchased  from  the  Danish  Hydraulic  Institute  and 
has  the  advantage  of  being  updated  regularly  with  new  processes  and  features.  In  this 
model  coastal  beds  can  be  rippled  or  nonrippled,  graded  or  of  uniform  material  and  surface 
waves  can  be  breaking  or  non-breaking  with  these  types  of  models  classified  as  being  the 
"best  way  forward  in  developing  more  advanced,  accurate,  reliable  sediment  transport 
prediction  methods"  (Soulsby,  1997). 

An  interesting  observation  of  these  types  of  models  was  made  by  Deigaard  (1998).  He 
observed  that  they  are  comparable  to  much  simpler  models  (namely  the  Bailard  formula). 
Hence  it  is  thought  that  the  more  complex  models  could  be  used  to  calibrate  the  simpler 
models,  which  have  the  major  advantge  of  not  needing  so  much  computing  resources. 

4.4  Summary 

From  critically  reading  the  available  literature  it  has  been  concluded  that  the  best  simple 
model  for  calculating  the  longshore  sediment  transport  rate  (or  the  vertically  integrated 
sediment  concentration  multiplied  by  longshore  velocity)  in  the  nearshore  (where  the  waves 
have  broken)  is  the  Queens  formula  (Holman,  1995)  with  the  model  by  Bailard  (1981)  used 
for  non-rippled  beaches  where  strong  external  currents  exist.  For  rippled  beds  the  models  by 
Soulsby  or  Van  Rijn  should  be  used.  Finally,  the  computer  intensive  STP  (Sediment 
Transport  at  a  Point)  model  by  tiie  Danish  Hydraulic  Institute  is  recommended  for 
calculating  the  sediment  transport  over  a  grid  and  can  be  applied  to  rippled  or  non-rippled 
beds  and  accounts  for  broken  as  well  as  xmbroken  waves.  This  model,  however,  is  expensive 
to  purchase. 
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5.  Vertical  Distribution  of  Suspended  Sediment 

The  vertical  distribution  of  sediment  concentration  in  the  water  column  within  and  seaward 
of  the  surf  zone  is  an  important  parameter  in  predicting  the  rmderwater  visibility.  By 
knowing  the  visibility  range,  well  informed  decisions  can  then  be  made  as  to  how  the 
detection  of  objects  within  the  nearshore  can  proceed. 

5.1  Suspended  Sediment  Concentrations  due  to  Currents: 

Due  to  the  large  tidal  amplitudes  which  occur  along  the  northern  Australian  coastline  (see 
Fig#22)  it  was  envisaged  that  there  may  be  occasions  when  strong  tidal  velocities  exist  in  the 
absence  of  any  significant  waves.  Hence,  it  was  thought  that  some  time  should  be  taken  to 
outline  those  models  which  calculate  the  vertical  suspended  sediment  concentration  arising 
solely  as  a  result  of  currents. 

In  the  sea,  the  model  which  appears  to  correspond  best  to  data  is  that  developed  by  Van  Rijn 
(Soulsby,  1997).  This  model  was  based  upon  the  steady  state  balance  between  the  settling  of 
sediment  due  to  gravity  and  the  upward  diffusion  resulting  from  turbulence.  This  balance 
can  be  expressed  by  equation#12,  with  Van  Rijn  making  the  assumption  that  the  turbulent 
diffusivity  of  sediment  varies  paraboUcaUy  in  the  lower  half  of  the  column  and  is  constant  in 
the  upper  half.  The  model  by  Van  Rijn  can  be  applied  to  regions  which  do,  or  do  not  have, 
rippled  beds. 


(O.C{z)=-K, 


dC(z) 

dz 


(12) 


where: 

1.  C(z)  is  the  suspended  sediment  concentration  at  a  height  z  above  the  bottom, 

2.  Ks  is  the  turbulent  diffusivity  of  sediment. 

Van  Rijn  Model: 


(13) 


With 


b  =  — +  B, 
B,  ' 
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„  ,  (0,  ^ 
B,  =1  +  2 

[u.) 

Bj  =2 

for  0.1  <1 

(13b) 

/  nO.8 

B^  =2.5  ^ 
[u.) 

Bj  =2 

/  f  Y  / 

for  0.01  <“/  <1 
^0.65  J  /u* 

for  ©5  >  u,  or  Zj  >  0.  IH 

(13c) 


where: 

1.  Za  is  a  reference  height  near  the  seabed, 

2.  Ca  sediment  reference  concentration  at  the  height  Za 

3.  H  is  the  water  depth 

4.  b  is  the  Rouse  number  ( — — )  or  suspension  parameter  (Soulsby,  1997) 

0.40u. 


5. 


is  the  friction  velocity 


)  where  is  the  vertical  stress  and  is  the  water-sediment 


mixture  density. 


The  reference  concentration  Ca  and  the  reference  height  Za  at  which  it  is  calculated  must  both 
be  evaluated  before  the  model  can  be  employed.  Several  formulae  have  been  developed 
with  Garcia  and  Parker  (1991)  testing  seven  and  recommending  the  following  two: 


Smith  and  McLean  (1977): 


0.00 156T, 
1-I-0.0024T, 


(13d)  at  a  height 


pg(s-l)  12 


(13e) 


where: 

(1)  = 

(2)  is  the  skin  friction  bed  shear  stress 

(3)  X„  is  the  threshold  bed  shear  stress  to  initiate  sediment  movement 
Van  Rijn  (1984) 

(14)  at  a  height 


Where 

(1)  A  5  is  the  height  of  sand  bed  waves  (see  page  116  of  Soulsby  (1996)  for  formulations) 
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(2)  D*  is  the  dimensional  grain  size  = 
method). 


(see  section  3.4  for  an  evaluation 


5.2  Suspended  Sediments  Concentrations  due  to  Unbroken  Waves: 

In  nearshore  regions  seaward  of  the  breaker  zone  the  suspended  sediment  resulting  from 
wave  induced  turbulence  is  largely  confined  to  the  relatively  thin  (few  cm's  or  mm's)  bottom 
botmdary  layer.  As  this  thickness  is  too  small  to  be  of  importance  to  divers  or  robotic 
vehicles,  this  section  has  been  included  for  purely  academic  reasons. 

For  rippled  seabeds  the  turbulent  sediment  diffusivity  is  assumed  constant  with  height  and 
the  concentration  profile  is  given  by  the  following  equation: 


(15) 


Where: 

1.  Co  is  a  reference  concentration  at  the  seabed  (z=0) 

2.  1  is  a  sediment  concentration  decay  length  scale, 

For  Co  and  1,  various  expressions  have  been  developed,  of  which  the  most  widely  accepted 
formulae  are  those  developed  by  Nielsen  (1992)  (Soulsby,  1997)  (see  page  145  of  Soulsby 
(1997)  for  a  complete  description). 

As  the  above  equation  was  developed  for  rippled  seabeds  (as  would  be  expected  in  the 
presence  of  surface  waves)  then  it  should  come  as  no  surprise  that  it  should  not  be  used  for 
sheet  flow  (flat-bed)  conditions.  In  such  situations  a  better  expression  to  use  is  analogues  to 
a  formula  which  is  sometimes  used  for  determining  suspended  sediment  concentrations  due 
to  currents  over  flat  beds.  The  formula  is  as  follows: 


(  ^ 

-b 

C(z)=C, 

z 

(16) 


Where: 

1.  b  is  a  constant  dependent  upon  the  median  grain  diameter. 

5.3  Suspended  Sediment  Concentration  due  to  Currents  and  Unbroken 
Waves 

The  suspended  sediment  concentration  due  to  waves  and  currents  first  originates  in  the 
wave  boundary  layer  and  then  diffuses  upward  through  the  water  column  as  a  result  of  the 
turbulence  of  the  current.  Both  of  these  processes  are  effected  by  the  interaction  of  the  wave 
boundary  layer  and  the  current  boimdary  layer,  with  the  concentration  profiles,  although  not 
being  well  established,  typically  described  for  flat  beds  by  the  following  formulae  (Soulsby, 
1997). 
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5.  Tn,gx  is  the  maximum  bed  shear  stress  in  one  wave  q^cle 

6.  is  the  mean  bed  shear  stress  in  one  wave  cycle 

7.  Ca  can  be  calculated  using  either  the  Van  Rijn  (1984)  formula  or  the  Smith  and  McLean 
(1977)  formula  outlined  above. 

Meanwhile,  if  the  concentration  profile  C(z)  is  required  above  a  rippled  seabed,  then  the  above 
model  can  stiU  be  appHed,  although  skin-friction  values  of  and  should  be  used  in  the  calculation 
of  the  reference  concentration  and  total  stress  values  should  be  used  to  calculate  bm  and  bmax . 

AU  of  the  above  formulations  for  the  suspended  sediment  profile  are  empirically  based  with 
each  based  upon  a  particular  set  of  assumptions.  The  most  advanced  method  for  calculating 
the  suspended  sediment  concentrations  (as  weU  as  transport  rates)  is  by  means  of  a 
numerical  model  of  the  wave  plus  current  botmdctry  layer,  utilising  some  form  of  turbulent 
energy  closure  scheme.  The  Danish  STP  model  is  one  such  model  of  this  type  and  was 
reported  above  for  longshore  sediment  transport.  It  has  an  advantage  that  it  can  be  applied 
to  plane  as  well  as  rippled  beds  (Soulsby,  1997). 

5.4  Suspended  Sediment  Concentration  under  Broken  Waves 

Waves  are  rarely  monochromatic  with  a  spectrum  of  wave  frequencies  and  heights  normally 
occurring.  As  waves  normally  break  along  any  shoreline  as  a  function  of  depth,  there  is  a 
horizontal  distance  between  where  the  first  waves  break  and  where  all  waves  have  broken. 
Such  a  region  is  called  the  "outer  zone"  of  the  surf  zone  (see  Fig#23)  (Fredsoe  and  Deigaard, 
1992),  wi&  the  complexity  of  the  turbulent  kinetic  energy  associated  with  breaking  of  waves 
in  this  region  not  sxifficiently  known  to  permit  the  suspended  sediment  concentrations  to  be 
predicted  (Fredsoe  and  Deigaard,  1992). 
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FigureU23  A  schematic  diagram  showing  how  the  surf  zone  can  be  disected  into  an  outer  zone  (where 
wave  breaking  occurs)  and  an  inner  zone  (where  all  broken  waves  have  organised 
themselves  as  hydraulic  bores).  Taken  from  Fredsoe  and  Deigaard  (1992). 

In  the  inner  zone,  however,  where  the  broken  waves  have  arranged  themselves  as  hydraulic 
bores,  the  effect  of  the  turbulence  upon  the  suspended  sediment  concentration  can  be 
described  by  using  the  turbulent  viscosity  formulation: 


(19) 


1.  1  is  the  length  scale  of  tiie  turbtilence  (see  page  166  of  Fredsoe  and  Deigaard  (1992)  for 
more  information  (and  equations)). 

2.  k  is  the  turbulent  kinetic  energy. 

The  turbulent  viscosity  formulation  can  then  be  placed  into  the  sediment  continuity  equation 
which  is  then  solved  numerically  to  obtain  the  suspended  sediment  concentration  (Fredsoe 
and  Deigaard,  1992). 


VT=lVk 

where: 


5.5  Summary 

For  forecasting  the  sediment  concentration  profile  in  the  vertical  (which  is  very  helpful  in 
forecasting  the  horizontal  visibility  range)  due  solely  to  the  action  of  horizontal  currents  the 
model  by  Van  Rijn  is  recommended  and  it  can  be  applied  to  regions  which  do,  or  do  not 
have,  rippled  beds.  The  sediment  concentration  in  the  vertical  due  to  the  action  of  currents 
and  unbroken  waves  is  easily  modelled  by  assuming  that  the  concentration  is  linearly  varied. 
The  sediment  concentration  due  to  broken  waves  can  only  be  forecasted  in  the  inner  surf 
zone.  This  is  because  it  is  only  here  that  the  broken  waves  have  arranged  themselves  as 
hydraulic  bores  with  the  effect  of  the  turbulence  upon  the  suspended  sediment  concentration 
described  by  the  hirbulent  viscosity  formulation  in  equation#19. 

5.6  Recommendations 

Although  numerical  models  are  available  which  predict  the  concentration  of  sediment  in  the 
vertical  under  a  wide  variety  of  environmental  conditions,  by  analysing  such  models  it 
should  quickly  be  seen  that  it  would  be  difficult  to  apply  them  in  an  operational  sense  by  the’ 
ADF.  The  reason  is  that  virtually  all  the  models  specified  above  require  a  lot  of  in-situ  data 
which  would  require  new  sktUs  to  be  learnt  by  ADF  personnel  and  considerable  man  hours 
expended  in  gathering  and  analysing  the  necessary  in-situ  data.  Hence,  none  of  the  above 
models  can  be  recommended  for  usage. 
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6.  Morphological  Modelling 

In  this  section,  models  which  allow  the  broad  scale  cross-shore  characteristics  of  the  beach 
morphology  to  be  predicted  are  examined  using  Short  (1999).  Also  examined  are  models 
that  predict  how  such  characteristics  vary  in  time.  These  temporal  models  have  been 
included  in  this  report  for  two  reasons.  First,  they  permit  the  determination  of  how  the 
topography  has  changed  due  to  storms.  Secondly,  they  permit  the  evolution  of  the  bottom 
topography  (due  to  changes  in  wave,  tide,  wind  and  current  conditions)  to  be  described  over 
longer  time  scales  (for  example  monthly  and  seasonal). 

6.1  Spatial  Indicative  Models 

From  reading  the  current  literature  it  was  concluded  that  there  are  currently  no  numerical 
models  in  existence  which  permit  the  two  dimensional  (longshore  and  cross-shore) 
morphology  of  the  surf  zone  to  be  predicted.  There  are  models  which  allow  the  user  to 
observe  what  morphology  will  result  when  different  combinations  of  edge  waves  exist  (for 
example  Holman  and  Bowen  (1982)),  but  such  models  are  usually  used  to  determine  which 
combinations  of  edge  wave  modes  must  be  present  in  the  surf  zone  to  form  a  particular  type 
of  morphology.  Such  models  cannot  be  used  in  a  predictive  or  operational  mode. 

There  are  models  however,  which  indicate  the  cross-shore  characteristics  of  morphology, 
such  as  the  number  of  longshore  bars,  where  the  bars  are  located  and  the  spacing  between 
any  rip  channels,  and  these  are  outlined  extensively  in  Short  (1999).  Also  outlined  in  Short 
(1999)  are  empirical  formulations  to  predict  the  slope  of  the  beachface. 

Predicting  which  beaches  will  have  longshore  bars  and  also  predicting  how  many  bars  will 
exist  parallel  to  any  beach  has  been  addressed  by  several  authors  with  Short  and  Aagaard 
(1993)  reviewing  them  and  using  field  data  from  microtidal  sandy  beaches  to  introduce  a 
new  bar  parameter  (BJ.  This  parameter  assumes  that  the: 

•  existence  of  any  longshore  bars  is  due  to  the  existence  of  a  dominant  cross-shore  standing 
infragravity  wave. 

•  cross  shore  profile  can  be  approximated  by  a  linear  slope  tan(p)  which  terminates  at  a 
constant  depth  at  a  distance  of  x^  from  tiie  shoreline. 

•  infragravity  wave  period  is  directly  proportional  to  the  incident  wave  period. 

The  equation  of  this  parameter  is  outlined  below  with  the  magnitude  of  it  versus  the  number 
of  bars  outlined  in  table  #1. 


B.  = 


Xs 

gT^  tan|3 


(20)  (Short  and  Aagaard,  1993). 


where: 

1.  T  is  the  incident  wave  period  and  used  as  a  surrogate  for  the  infragravity  wave  period. 
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Tableitl  The  number  of  bars  which  correspond  to  the  magnitudes  ofB^.  Taken  from  Short  (1999). 


<20 

20-50 

50-100 

100-400 

>400 

Number  of  bars 

0 

1 

2 

3 

4 

For  longshore  bars  which  have  formed  due  to  the  presence  of  standing  infragravity  waves 
the  location  of  the  sand  bars  crests  or  troughs  in  the  offshore  direction  can  be  predicted.  This 
is  because  such  features  will  either  be  at  the  standing  infragravity  wave  nodes  or  antinodes 
when  the  mode  of  sediment  movement  is  by  bedload  or  suspended  load,  respectively.  The 
distance  (x)  from  the  shoreline  to  consecutive  nodes  of  the  standing  wave  is  given  by; 


gTe  tanpx 
471^ 


(21)  (Bowen,  1980) 


where: 

1.  Te  is  the  period  of  the  standing  infragravity  wave 

2.  X  is  the  non-dimensional  distance  from  the  shoreline  to  a  given  node  or  antinode  and  is 
given  in  Table#2. 

Table#!  Non-dimensional  distance  from  the  shoreline  to  nodes  of  infravravihi  waves.  Taken  from 
Short  (1999). 


Wave  Type 

Node  1 

Node  2 

Nodes 

Node  4 

Node  5 

Edge  Wave  Model 

1.5 

Edge  Wave  Mode2 

1.5 

8.5 

Edge  Wave  Mode3 

1.5 

8.0 

22.0 

Edge  Wave  Mode4 

1.5 

7.9 

20.4 

42.3 

Edge  Wave  Mode5 

1.4 

7.8 

19.8 

39.0 

69.5 

Leaky  Wave 

1.4 

6.6 

17.6 

33.8 

54.8 

Antinode 

Antinode 

Antinode 

Antinode 

Antinode 

1 

2 

3 

4 

5 

Edge  Wave  Model 

4.5 

Edge  Wave  Mode2 

3.9 

16.1 

Edge  Wave  Mode3 

3.8 

13.6 

35.1 

Edge  Wave  Mode4 

3.7 

13.0 

29.6 

61.6 

Edge  Wave  Mode5 

3.7 

12.8 

28.1 

52.2 

95.7 

Leaky  Wave 

3.5 

12.3 

25.8 

44.5 

67.8 

For  standing  edge  waves  in  the  longshore  direction  the  continuation  of  any  longshore  bars 
wiU  be  terminated  by  the  formation  of  rip  channels.  The  location  of  such  channels,  although 
being  at  the  nodes  of  the  longshore  standing  edge  wave,  cannot  be  predicted.  However,  the 
distance  between  them  can  be  calculated  as  it  is  half  the  wavelength  of  the  edge  waves  and 
can  be  obtained  from  the  following  formula. 


^^AT2(2n  +  l)tanp 
2  47t  ^  ^ 

where: 


(22)  (Short,  1999) 
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1.  is  the  edge  wave  length. 

2.  n  is  the  mode  number  of  the  edge  wave. 

This  equation  however,  can  be  simplified  for  steep  reflective  beaches  (tanp  >0.01).  This  is 
because  for  such  beaches,  Guza  and  Davis  (1974)  made  the  observation  that  mode  0 
subharmonic  edge  waves  (edge  waves  having  half  the  period  of  the  incident  waves) 
predominate  with  smaller  synchronous  edge  waves  (edge  waves  having  the  same  period  as 
the  incident  waves)  also  having  a  tendency  to  dominate  on  such  beaches.  Using  these 
observations  the  distance  between  rip  channels  dissecting  any  longshore  bar  can  be 
simplified  to  the  following  equations: 


2  k 


(23a)  subharmonic  edge  waves 


^  =  -i-T^2^anp 

2  2k  ' 


(23b)  synchronous  edge  waves 


In  the  above  paragraphs  regarding  predicting  which  beaches  will  have  bars  formed  by  cross¬ 
shore  standing  infragravity  waves,  how  many  will  form  and  their  associated  offshore  crest 
positions,  it  is  assumed  in  the  models  that  there  is  a  dominant  infragravity  frequency.  This 
assumption  however  is  rarely  vahd  with  the  infragravity  energy  spectrum  being  broad 
banded  and  often  with  no  dominant  modes.  Only  in  storm  conditions  when  the  surf  zones 
are  wider  has  dominant  infragravity  modes  been  observed  in  the  surf  zone  (Short  and 
Aagaard,  1991).  Even  with  this  assumption  not  being  valid  most  of  the  time  however.  Short 
and  Aagaard  found  that  the  ability  of  to  predict  the  number  of  likely  bars  was  siuprisingly 
good  (Ibid). 

For  single  longshore  bars  which  have  formed  as  a  result  of  the  convergence  of  the  onshore 
wave  induced  movement  of  sand  with  the  offshore  movement  caused  by  pltmging  waves 
breaking  Holman  and  SaUenger  (1993)  have  constructed  a  model  which  predicts  the  location 
of  the  bar  crest  at  Xbar  meters  from  the  shoreline  with  Xbar  taking  tire  following  form: 


_  hbar  _  yHb 
tanP  tanP 


(24) 


where: 

1 .  hbar  is  the  depth  where  the  waves  break. 

2.  Hb  is  the  wave  breaker  height. 

3.  Yis  a  constant  of  order  1  and  expresses  the  empirical  ratio  of  wave  height  to  depth  of 
breaking  waves.  Thus  wave  height  provides  the  primary  scaling  for  bar  distance. 


One  drawback  of  the  above  model  is  that  waves  are  rarely  monochromatic  but  instead  have 
a  spectrum  of  wave  heights.  Thus  instead  of  monochromatic  waves  breaking  at  a  single 
point  and  generating  a  sand  bar,  normal  waves  break  over  a  region  and  it  is  unclear  if  a  bar 
will  form  by  this  mechanism  or  if  one  will  form  where  it  will  be  located.  However,  it  is 
known  that  once  a  bar  forms,  it  does  have  the  ability  to  focus  wave  energy  which  allows  it  to 
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grow  further  through  positive  feedback  (Holman  and  Sallenger,  1993).  Similar  positive 
feedback  mechanisms  have  been  observed  between  pre-existing  multiple  bars  and  cross¬ 
shore  standing  edge  waves  (Ibid). 

The  slope  of  the  beachface  or  swash  zone  is  determined  from  the  balance  between  two 
factors;  the  strength  of  the  swash,  which  results  in  onshore  sediment  movement,  and  gravity 
which  results  in  offshore  sediment  movement.  The  region  over  which  the  beachface  slope  is 
defined  is  within  the  mid-intertidal  zone  with  the  beachface  slope  normally  being  concave  in 
the  seaward  direction  and  rarely  being  planar.  For  a  net  offshore  sediment  movement  the 
beachface  slope  becomes  flatter,  whilst  the  beachface  becomes  steeper  when  there  is  a  net 
onshore  movement  of  sediment. 

Sunamura  (1984b)  using  the  results  from  over  35  publications  in  conjunction  with 
dimensional  analysis  derived  the  following  empirical  formula  for  the  equilibrium  beach  face 
slope  in  unsaturated  conditions. 


0.12 

(25) 


where: 

1.  P  is  the  beach  face  slope  in  degrees. 


In  obtaining  this  relationship  it  was  observed  that  there  was  a  lot  of  scatter  of  the  data  about 
this  function.  Figure#  24  shows  this  function  and  the  data  it  was  obtained  from.  From  this 
observation  it  was  concluded  that  there  were  other  factors  not  considered  in  equation  #25 
which  have  an  influence  of  the  slope  of  the  beach  face. 


Figure#24  Data  points  generated  when  the  measured  beachface  slope  is  plotted  against  the 


dimensionless  parameter 


The  least  square  fit  is  also  shown.  From  Short  (1999). 


One  parameter  which  can  have  an  influence  on  the  slope  of  the  beach  face  is  the  level  of  the 
water  table  at  the  beachface  or  the  water  content  of  the  beachface  sediment.  Such  a  factor  is 
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partioilarly  important  on  macro-tidal  beaches  where  two  different  intertidal  profiles  can 
exist.  In  a  effort  to  extend  the  work  of  Sunamura  to  include  saturated  beaches.  Turner  (1994) 
investigated  15  macro-tidal  beaches  m  Queensland,  Australia  where  there  were  lower  parts 
of  the  beachface  that  were  permanently  saturated.  By  re-analysing  the  work  of  Suncunura, 
Kriebel,  Kraus  and  Larson  (1991)  derived  the  following  relationship  (equation#26)  for  the 
eqiiilibrium  slope  of  saturated  beachfaces.  It  should  be  noted  however,  Aat,  in  comparison 
to  the  number  of  data  points  used  to  derive  the  imsaturated  beachface  slope,  only  a  limited 
munber  were  used  to  derive  the  following  equation. 


p  =  0.05 


(26) 


where: 

all  parameters  have  been  previously  defined. 

In  order  to  utilise  equation#26  the  settling  velocity  wotdd  need  to  be  either  measured  or 
calculated.  Section#3.4  outlines  the  way  the  ADF  could  calculate  this  value. 

6.1.1  Summary 

Currently  it  is  possible  to  predict  the  number  of  longshore  bars  that  exist  along  any  beach,  the 
location  of  the  crests  in  the  offshore  direction,  and  the  beachface  slope  for  saturated  and 
unsaturated  beachfaces.  It  is  not  possible  to  use  munerical  models  to  predict  the  longshore  and 
cross-shore  two  dimensional  morphology  within  the  surf  zone  at  any  point  in  time.  Such 
information  must  be  collected  using  instruments  such  as  Laser  Airborne  Depth  Sounders. 

6.1.2  Recommendations 

From  analysing  the  constituents  of  the  above  models  and  comparing  them  to  what  the  ADF 
is  likely  to  be  able  to  measure  or  predict,  it  is  foreseen  that  it  would  be  possible  for  the  ADF 
to  use  the  Bar  Parameter  (B*)  to  predict  which  beaches  wiU  have  bars  forming  by  cross-shore 
standing  infragravity  waves  and  approximately  how  many  are  likely  to  form.  The  offshore 
position  of  the  bar  crests  formed  by  infragravity  waves  as  well  as  the  longshore  rip  channels 
formed  by  the  presence  of  edge  waves,  however,  cannot  be  predicted  as  the  dominant 
infragravity  wave  period  needs  to  be  known  and  such  a  parcuneter  cannot  be  predicted  and 
is  very  difficult  to  measure. 

For  single  longshore  bars  which  have  formed  by  monochromatic  breaking  waves,  the  model 
which  describes  the  offshore  location  of  the  bars  crests  (equation#24)  should  be  viable  to  use 
by  the  ADF,  as  it  is  only  requires  very  basic  input  parameters. 

The  models  described  above  for  predicting  the  slope  of  the  beachface  imder  saturated  and  - 
non-saturated  conditions,  can  be  used  by  the  ADF,  although  the  model  for  saturated 
beachfaces  requires  that  a  small  sediment  sample  be  obtained  and  the  median  grain  size 
visually  estimated. 
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6.2  Temporal  Evolution  Models 

6.2.1  Introduction 

Broadly  speaking,  there  are  three  types  of  morphological  models  which  predict  how 
different  aspects  of  the  surf  zone  morphology  change  in  time.  These  include  beach 
planshape  models,  coastal  profile  models,  and  coastal  area  models  (Soulsby,  1997).  The 
beach  planshape  models  predict  the  evolution  of  the  shape  of  the  shoreline  (for  example  the 
still  water  level).  They  achieve  this  by  first  calculating  longshore  transports  at  equally 
spaced  points  along  a  shore  and  by  then  applying  the  sediment  budget.  These  types  of 
models,  however,  work  over  time  scales  of  years  and  decades.  Hence,  they  are  not  relevant 
from  the  point  of  view  of  nearshore  modelling  for  supporting  amphibious  operations.  The 
other  two  types  of  models  are  relevant  to  the  support  of  amphibious  operations  and  are 
explained  in  the  following  paragraphs. 

6.2.2  Coastal  Profile  Models 


Coastal  profile  models  are  only  concerned  with  how  the  cross-shore  beach  profile  (or 
topography)  evolves  in  time  due  to  varying  environmental  forces.  These  types  of  models 
are  more  computationally  demanding  than  beach  planshape  models  and  as  such  can  only  be 
apphed  to  smaller  time  scales  (typically  several  weeks  or  months).  This  type  of  modelling 
can  be  achieved  by  numerous  methods  and  include:  (1)  Descriptive  models,  (2)  Equilibrium 
models,  (3)  Empirical  models,  and  (4)  Process  based  models,  (Roelvink  and  Broker  (1993). 


Descriptive  models  are  useful  in  predicting  the  typical  topography  expected  at  any  beach  as 
weU  as  transitions  between  different  beach  states.  Such  predictions  are  achieved  by  using 
governing  parameters  as  indicators,  with  such  parameters  determined  from  beach 
observations  obtained  over  as  long  a  period  as  possible  (preferably  several  months  or  more). 
A  common  governing  parameter  is  Dean's  (1973)  parameter  Q(Dean,  1991;  Roelvink  and 
Broker,  1993;  Short,  1999),  which  was  called  the  dimensionless  fall  velocity  in  section  3.4  (see 
Eqn.#9).  This  parameter  is  based  on  the  observation  that  wave  dominated  beaches  can  be 
classified  using  wave  height  (Hb),  wave  period  (T)  and  grain  size  (as  defined  by  the  sediment 
fall  velocity  CQj )  (Short,  1999).  Descriptive  models,  however,  cannot  be  used  to  infer  beach 
changes  caused  by  man  made  (or  proposed  man-made)  structures.  They  also  have  limited 
quantitative  capabilities  (Roelvink  and  Broker,  1993). 


Dean's  parameter 


Equilibrium  profile  models  are  basic  algebraical  models  which  predict  the  new  cross-shore  beach 
profile  that  would  occur  if  aU  environmental  forces  were  in  equilibrium  (which  would  be  when 
rate  of  sedimentation  equals  the  rate  of  erosion).  Equilibrium  profile  models  are  useful  in  cases 
of  sandy  beaches  where  longshore  sand  transport  rates  are  negligible  and  can  be  used  for 
establishing  equilibrium  profiles  at  time  scales  from  those  of  long  storms  to  years. 


One  of  the  most  widely  used  equilibrium  profile  models  is  that  by  Dean  (1977)  (as  cited  by 
Roelvink  and  Broker,  1993).  This  model  describes  the  equilibrium  cross-shore  profile  for 
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monotonic  profiles  that  would  occur  because  of  a  rise  in  sea-level  due  to  such  influences  as 
storms.  The  model  is  a  2/3  power  shape  law  with  a  constant  of  proportionality  (A)  that  is 
dependent  upon  the  sediment  distribution  (see  Eqn.#28).  Like  equilibrium  models  in 
general,  it  assumes  ttiat  there  is  no  gradient  in  the  longshore  sand  transport  Unfortunately, 
this  model  does  not  accoimt  for  variations  in  the  wave  conditions,  wind  speeds  or  currents 
associated  with  varying  sea  levels. 


Dean's  Rule 


z=Ax' 


(28) 


Where 

1.  z  and  X  represent  the  vertical  and  offshore  axis,  respectively. 

2.  A  is  a  scale  factor  which  is  primarily  dependent  upon  sediment  characteristics. 

Empirical  profile  evolution  models  illustrate  tiie  evolution  of  the  cross-shore  beach 
topography  towards  its  equilibrium  profile  and  are  basic  numerical  models  in  the  sense  that 
they  have  empirically  derived  constants  and  relationships.  These  models  can  be  applied  to 
monotonic  beach  profiles  as  well  as  beaches  with  longshore  bars  and  troughs,  and  have  an 
important  application  of  evaluating  damage  caused  by  storms  of  limited  duration  (Roelvink 
and  Broker,  1993).  One  of  the  disadvantages  of  these  types  of  models  is  that  empirical 
coefficients  must  be  determined  for  each  beach. 


Process-based  models  are  models  which  have  a  module  for  each  process  which  contributes 
to  the  morphological  changes.  The  modules  can  be  for  tides,  waves,  nearshore  velocities, 
cross-shore  sediment  transport  and  morphology,  with  the  morphological  modules  normally 
based  upon  solving  the  sediment  budget  equation  for  bed  height  (zb).  The  sediment  budget 
equation  for  non-horizontal  beds  is  given  by  the  following  equation  (taken  from  de  Vriend 
ef-flZ.,  1993). 
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where: 

1.  t  is  time  (sec), 

2.  Ep  is  the  bed  porosity  (dimensionless), 

3.  s  and  n  are  tire  distance  along  and  orthogonal  to  streamlines  (m), 

4.  qj  is  the  total  sediment  transport  for  a  horizontal  bed, 

5.  bi  and  hi  are  constants, 

6.  h  is  the  water  depth  (m), 

7.  Rn  is  the  radius  of  curvature  of  the  streamlines, 

8.  Is  and  L  are  metric  coefficients  related  to  the  curvilinear  coordinate  system, 

9.  P  is  the  downslope  (gravitational)  sediment  transport. 

Process  based  models  are  computer  extensive  models  and  have  the  major  disadvantage  that 
they  require  knowledge  of  the  nearshore  sediment  characteristics  to  be  applied.  They  also 
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have  the  disadvantage  that  they  should  only  be  applied  by  personnel  with  a  strong 
background  in  nearshore  modeUing. 

6.2.3  Coastal  Area  Models 

These  two  dimensional  (area)  morphology  models  are  essentially  the  same  as  the  process- 
based  models,  with  the  only  difference  being  that  the  longshore  sediment  transport  as  well 
as  the  cross-shore  transport  are  included  in  solving  the  sediment  budget  equation  for  bed 
depth.  These  types  of  models  have  individual  modules  (or  numerical  models)  for  waves, 
currents  and  sediments  and  in  the  long  term  are  thought  to  be  the  best  way  to  model  aU 
aspects  of  the  nearshore  environment.  However,  they  are  computer  resource  expensive 
models,  can  cost  a  large  amoimt  of  money  (>  $100K  AUS)  and  need  to  be  adjusted  for  each 
apphcation  by  experienced  personnel.  They  are  normally  purchased  as  whole  packages  with 
associated  software  from  such  organisations  as  HR  Wallingford,  Wallingford,  UK;  cind  Delft 
Hydraulics,  Netherlands. 

6.2.4  Summary 

From  reading  the  hterature  it  has  been  observed  that  there  are  three  types  of  morphological 
models  which  predict  how  the  morphology  of  the  surf  changes  in  time.  These  include,  beach 
planshape  models,  coastal  profile  models  and  coastal  area  models  (Soulsby,  1997).  Although 
coastal  profile  models  and  coastal  area  models  were  deemed  useful  in  the  planning  of 
amphibious  operations,  beach  planshape  models,  were  not.  This  is  because  such  models 
predict  the  evolution  of  the  shape  of  the  shoreline  over  time  scales  of  years  and  decades. 

6.2.5  Recommendations 

By  comparing  the  complexity  of  the  temporal  morphology  models  described  above  against 
the  in-situ  parameters  the  ADF  would  need  to  measure  to  use  them,  it  was  quickly  concluded 
that  the  ADF  was  not  in  a  position  to  use  such  models. 

7.  Summary 

Before  summarising  the  models  foimd  within  the  scientific  literature  that  were  deemed 
viable  for  the  ADF,  it  should  be  stated  that  in  numerous  cases  the  best  predictive  models 
could  not  be  recommended  because  they  were  either  too  complicated  or  required  input 
parameters  which  were  too  difficult  to  acquire.  Hence  a  subliminal  lesson  from  the  work 
conducted  here,  is  that  there  can  sometimes  be  a  large  difference  between  what  is  available 
from  the  scientific  community  as  a  result  of  research  and  what  can  be  viably  implemented 
from  an  operational  perspective. 

Keeping  the  above  paragraph  in  mind,  it  was  concluded  that  the  most  useful  model  for 
calculating  the  longshore  velocity  within  the  surf  zone  as  a  result  of  breaking  waves  was  that 
developed  by  Komar  and  Inman  (1970).  Such  a  model  is  only  applicable  to  monotonically 
sloping  beaches  of  constant  slope  and  for  monochromatic  waves,  is  not  as  accurate  as  more 
complicated  models  although  was  recommended  because  it  has  the  large  benefit  of  only 
requiring  wave  characteristics  as  input. 
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Currently  there  are  no  numerical  models  in  existence  which  permit  the  longshore  and 
offshore  currents  in  the  surf  zone,  generated  by  gradients  in  radiation  stress,  to  be  predicted 
as  a  function  of  time  and  space.  Such  models  would  require  a  large  number  of  parameters  as 
input  and  it  is  thought  that  such  models  will  not  be  available  for  at  least  5  years.  What  has 
been  foxmd  are  the  range  in  velocity  strengths  of  rips  and  mean  values  in  longshore  spacing 
of  rips.  A  quantitative  model  developed  by  Martens  et-al  was  also  found  which  permits  the 
ability  of  topographic  rips  to  form  in  embayments  to  be  predicted. 

Models  for  forecasting  of  the  sediment  concentration  m  the  vertical  (which  is  very  helpful  in 
forecasting  the  horizontal  visibility  range)  due  to  the  action  of  currents,  tmbroken  and 
broken  waves  were  outlined  in  the  text.  However,  it  was  quickly  concluded  that  due  to  such 
models  requiring  complex  parameters  as  input  none  could  be  recommended  to  the  ADF. 

For  predicting  nearshore  morphology  characteristics  it  was  found  that  there  are  no 
numerical  models  in  existence  which  permit  such  characteristics  to  be  predicted  in  the 
longshore  and  offshore  direction.  Such  information  would  need  to  be  gathered  from  high 
resolution  instruments  as  Laser  Airborne  Depth  Sounders.  What  was  foimd  were  models 
that  allowed  the  characteristics  of  longshore  bars  to  be  predicted.  In  particular  it  was 
concluded  that  it  would  be  possible  for  the  ADF  to  use  the  Bar  Parameter  (B.)  to  predict 
which  beaches  will  have  bars  forming  by  cross-shore  standing  infragravity  waves  and 
approximately  how  many  are  hkely  to  form.  Models  describing  the  offshore  position  of  the 
bar  crests  formed  by  standing  infragravity  waves  as  well  as  the  separation  distance  of 
longshore  rip  channels  were  also  identified,  however,  they  could  not  be  recommended  to  the 
ADF  as  the  dominant  infragravity  wave  period  needed  by  them  cannot  be  predicted  and  is 
difficult  to  measure. 

For  single  longshore  bars  which  have  formed  by  breaking  waves,  the  model  which  describes 
the  offshore  location  of  the  bars  crests  (equation#24)  should  be  viable  to  use  as  it  is  only 
requires  basic  input  parameters.  Lastly,  models  enabling  the  slope  of  the  beachface  to  be 
predicted  for  saturated  and  non-saturated  conditions  were  also  foimd  (equations#25  and#26) 
and  they  too  can  be  recommended  to  the  ADF,  but  only  if  a  small  sediment  sample  is 
obtained  and  the  median  grain  size  estimated. 

Finally  numerical  models  were  identified  within  the  literature  which  enable  the  temporal 
changes  that  occur  to  nearshore  morphology  from  time  scales  of  days  to  years  to  be 
predicted.  Such  models  however  were  foreseen  as  not  viably  useful  to  the  ADF  because  they 
require  in-situ  data  which  would  be  difficult  to  obtain. 

8.  Summary  Table 

Longshore  velocity  strength  on  beaches  with  no  rip  cells. 


Reference 

Assumptions 

Input  Parameters 

Komar  and  Inman  (1970) 

Plane  beaches  with 
monochromatic  waves 

Breaking  wave 
height,  wave 
propagating  angles 
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Thornton  and  Guza  (1986) 

Plane  beaches  with  irregular 
waves,  wave  length  to  depth 
ratio  <1/20 

Wave  propagating 
angles,  wave  type, 
water  depth,  surf 
zone  slope,  deep 
water  wave  height, 
peak  wave 
frequency,  coastline 
orientation. 

Smith,  Larson  and  Krauss 
(1993) 

Non-monotonic  beaches  with 
monochromatic  waves. 

Breaking  wave 
height,  propagation 
angle,  wave  period. 

Longshore  velocity  with  offshore  velocity  (two  dimensional  models) 


Prognostic  Parameter 

Reference 

Input  Parameters 

Longshore  spacing  between  edge 
wave  induced  rips 

Short  (1999) 

Edge  wave  period,  edge  wave 
node,  surf  zone  slope 

Longshore  velocity  variation  in 
the  longshore  direction  in  the 
trough  separating  a  bar  (dissected 
by  rip  channels)  and  the  swash 
zone. 

Dalrymple 

(1978) 

Cross-sectional  area  of 
trough,  offshore  wave 
steepness,  depth  over  bar 

Mega  Rip  formation 

Martens, 
Williams  and 
CoweU  (1999) 

Embayment  Dimensional 
Characteristics,  breaking 
wave  height,  wave  period 
and  sediment  settUng 
velocity. 

Vertical  distribution  of  suspended  sediment. 


Forcing  Parameter 

Reference 

Input  Parameters 

Horizontal  currents 

Van  Rijn 

Water  depth,  water  density, 
sediment  reference 
concentration,  vertical  stress, 
sediment  settUng  velocity. 

Horizontal  currents  with 
unbroken  waves 

Souslby  (1997) 

Mostly  the  same  as  the  above 
parameters. 

Broken  waves  in  the  inner  surf 

zone 

Fredsoe  and 

Deigaard 

(1992) 

Mostly  the  same  as  the  above 
parameters. 

Morphology  Characteristics  (spatially  indicative  models) 


Prognostic  Parameter 

Reference 

Input  Parameters 

Offshore  location  of  single  sand 
bars  formed  by  the  action  of 
breaking  waves. 

HoUman  and 
Sallenger  (1993) 

Wave  height  at  breaking  and 
bed  slope. 
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Bar  Parameter  (B*)  [if  bars  will 
form  by  reflected  infragravity 
waves,  and  if  so  how  many]. 

(Short  and 
Aagaard  (1993) 

Incident  wave  period,  surf 
zone  bed  slope,  distance  from 
shoreline  where  monotonic 
slope  terminates. 

Offshore  location  of  sand  bars 
formed  due  to  standing 
infragravity  waves 

Bowen  (1980) 

Infragravity  wave  period,  surf 
zone  bed  slope. 

Longshore  Spacing  between  rip 
channels. 

Short  (1999) 

Edge  wave  period,  edge  wave 
node  and  surf  zone  slope. 

Beach  face  slope 

Sunamura 

(1984b) 

Wave  height  at  breaking,  wave 
period  and  median  grain  size. 

Morphology  Qmaracteristics  (Temporally  indicative  Models) 


Type  of  Model 

Reference 

Input  Parameters 

Coastal  Profile  Models 

Dean  (1977) 

Breaking  wave  height, 
wave  period  and  sediment 
settling  velocity 

Coastal  Area  Models 

Incident  wave 
characteristics,  surf  zone 
sediment  characteristics, 
surf  zone  bathymetry  and 
tidal  characteristics. 
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